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ABSTRACT 

We present the discovery of a kinematically-cold stellar population along the southeastern minor axis 
of the Andromeda galaxy (M31) that is likely the forward continuation of M31's giant southern stream. 
This discovery was made in the course of an on-going spectroscopic survey of red giant branch (RGB) 
stars in M31 using the DEIMOS instrument on the Keck II 10-m telescop^EJ Stellar kinematics are 
investigated in eight fields located 9-30 kpc from M31's center (in projection). A likelihood method 
based on photometric and spectroscopic diagnostics is used to isolate confirmed M31 RGB stars from 
foreground Milky Way dwarf stars: for the first time, this is done without using radial velocity as 
a selection criterion, allowing an unbiased study of M31's stellar kinematics. The radial velocity 
distribution of the 1013 M31 RGB stars shows evidence for the presence of two components. The 
broad (hot) component has a velocity dispersion of tXy ph =129 km s -1 and presumably represents M31's 
virializcd spheroid. A significant fraction (19%) of the population is in a narrow (cold) component 
centered near M31's systemic velocity with a velocity dispersion that decreases with increasing radial 
distance, from c™ b= 55.5 km s^" 1 at i? pro j = 12 kpc to <J® = 10.6 km s" 1 (an intrinsic velocity 
dispersion of 9.5 km s _1 after accounting for velocity measurement error) at -R pro j = 18 kpc. The 
spatial and velocity distri bution of the cold component matches that of the "Southeast shelf" predicted 
by the lFardal et al.l (|2007l ) orbital model of the progenitor of the giant southern stream. The metallicity 
distribution of the cold component matches that of the giant southern stream, but is about 0.2 dex 
more metal rich on average than that of the hot spheroidal component. We discuss the implications 
of our discovery on the interpretation of the intermediate-age spheroid population found in this region 
in recent ultra-deep HST imaging studies. 

Subject headings: galaxies: substructure — galaxies: halo — galaxies: individual (M31) — stars: 
kinematics — techniques: spectroscopic 



1. INTRODUCTION 



Electronic address: kgilbert@ucolick.org, 

fardal@fcraol.astro.umass.ee u, jkalirai@ucolick.org, 
raja@ucolick.org, marla.geha@nrc-cnic.gc.ca, jcisler@ucolick.org, 
srm4n@virginia.edu, jo stheim@alumni.virginia.edu, 

rjpOi@virginia.edu, reitzel@astro.ucla.edu, ekirby@ucolick.org, 
cooper@astron.berkeley.edu 

1 UCO/Lick Observatory, Department of Astronomy & Astro- 
physics, University of California Santa Cruz, 1156 High Street, 
Santa Cruz, California 95064. 

2 Department of Astronomy, University of Massachusetts, 
Amherst, Massachusetts 01003. 

3 Hubble Fellow. 

4 NRC Herzberg Institute of Astrophysics, 5701 West Saanich 
Road, Victoria, British Columbia, Canada V9E 2E7. 

5 Plaskett Fellow. 

6 Fisk University /Vanderbilt University, Nashville, Tennessee 
37325. 



7 Visiting Student, University of California Santa Cruz. 

8 Department of Astronomy, University of Virginia, 



PO 



Box 400325, Charlottesville, VA 22904-4325. 

9 Department of Physics &; Astronomy, Knudsen Hall, Univer- 
sity of California, Los Angeles, California 90095. 

Department of Astronomy, Campbell Hall, University of 
California, Berkeley, California 94720. 

11 Data presented herein were obtained at the W. M. Keck Ob- 
servatory, which is operated as a scientific partnership among the 
California Institute of Technology, the University of California and 
the National Aeronautics and Space Administration. The Obser- 
vatory was made possible by the generous financial support of the 
W. M. Keck Foundation. 



In the current paradigm of hierarchical galaxy forma- 
tion, massive galaxies are built u p through a series of 
major and minor m erger events (|Searle fe Zinnl 119781 : 
I White fc Reesl fl978l ) . Observations of galaxies at high 
redshift show that merging systems are c ommon (e.g., 
lAbraham et al.lll996l: iLeFevre et alj|2000t Conselice et 
al. 120031 : ILo~ et al J 120071 ). and large-scale simulations 
of galaxy formation in a cosmological context have suc- 
cessfully reproduced many of the o bserved properties of 
galaxies and galax y clusters (e.g.. ISpringel et ail 120051 : 
iCroton et aLll2006l ). 

A consequence of hierarchical galaxy formation is that 
galactic stellar halos should be at least partially com- 
posed of the tidal debris from past accretion events. Nu- 
merical simulations and semi-analytic models of stellar 
halo formation have made great strides toward under- 
standing the properties of halos that are built up through 
tidal stripping of mer g ing systems (e.g. , I Johnston et al] 
[19961 : iJohnstod fl998l: iHelmi fc White) 119991: Helmi fc 
de Zeeuw l2000l : lBullock et al.ll2001l : lBullock fe Johnston! 
[III). Detailed comparisons between observations and 
simulations are needed to determine the fraction of stel- 
lar halos that are made up of tidal debris and to better 
understand the formation of galaxies in general. 

Recent discoveries of tidal streams in the stellar halos 
of the Milky Way (MW) and Andromeda (M31) galaxies 
are providing the most direct and detailed observational 
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constraints on theories of stellar halo formation. Among 
the most prom inent of these subs t ructures are the Sagit- 
tarius stream (llbata et alj 11994 iMaiewski et al.ll2003l ; 
iNewberg et all l2003f) . the Monoceros stream ( Yannv 
et al. 120031 : iRocha-Pinto et all 1200311 and the Maeel- 
lanic stream (Mathcws on et all 11974) in the MW, a nd 
the giant southern stream (GSS: llbata eTaDdOOiB) m 
M31. Additional substructures have been identified in 
M31 that are also likel y remnants of past mergers, such 
as th e Northeast shelf ([ Ferguson et al. 2002; Ibat a et all 
12004 iFardal et al.ll2006f ). a secondary cold component in 



the same physical lo cation as the giant southern stream 
(jKalirai et al.ll2006al ). and the various substructures iden- 
tified with the disk of M31 (|Ibata et al.ll2005l l Tidal dis- 
ruption has also bee n observed in M3 1's satellite galaxies 
M32 and NGC 205 (|Choi et al.H2002T ). 

In addition to providing insight into theoretical mod- 
els of stellar halo formation, the observed properties of 
tidal streams can be used to constrain the the galac- 
tic potential in which they are fou nd if sufficient phase- 
space information is available (e.g., Uohnston et al.lll999l . 
|2002| ; [Penarrubia et al.1 12006T ). Several attempts have 
been made to model the mass distribution of the MW 
using observed substructure, most of which have focused 
on the orbital properties of the Sagittarius stream (e.g., 
Ibata et al1l2001al : fHclmi 2004; Martuicz -Delgado et all 



bata ct al. 2004: Guhathakurta ct al. 2006: Kalirai 



2006al l The availability of this phase-space infor- 



200 4 Uohnston et al.ll2005l: lLaw et "all 12005k Fellhauer 
et al. 120061 1 

The GSS has been the focus of detailed modeling in 
M31, spurred on by recent observations. Imaging and 
photometry have revealed the physical extent o f the GSS 
([Ferguson et al.ll2002l:lMcConnachie et al.ll2003l: Fer guson 
et al. I2006D and pr ovided line-of-sight dista nces at vari- 
ous points along it (|McConnachie et aTll2003l ). while spec- 
troscopy has yielded the mean line-of-sight velocity and 
veloc ity dispersion of stream stars as a function of posi- 
tion { 
et al. 

mation has motivated sever al groups to model the orbit 
of the progenitor of the GSS (|Ibata et aI1 l2004: Fo nt et all 
120061: IFardal et all 120061 . |2007|) However, IFardal et all 
(2006) concluded that the degree to which the GSS can 
be used to constrain M31's mass distribution is limited 
by the current measurement uncertainties in the distance 
to the stream and the lack of a clearly identified compact 
stellar concentration that might correspond to the dense 
remnant core of the stream's progenitor galaxy. Further 
observational constraints on the orbit of the progenitor 
of the GSS, such as the identification of tidal debris from 
other pericentric pas sages, are needed to make progress. 

Towards this end, IFardal et all (|2007l . hereafter F07) 
show that several of the observed features in M31 can 
be explained as the forward continuation of the GSS. 
Their model makes predictions that can be tested by 
observations, including the stellar velocity distributions 
in the Northeast and Western shelves and the presence 
of a weaker shelf on the eastern side of the galaxy. This 
last shelf is expected to be most easily visible near the 
southeastern minor axis of M31 12 . 

12 Although this shell feature is expected to span an ~ 180° 
range in position angle, covering the eastern half of the galaxy, 
it is expected to be most easily observable in the southeast due 
to overlap with th e m uch denser Northeast shelf and M31's disk 
elsewhere (see Fig. [12} ■ 



This paper presents new substructure along M31's 
southeastern minor axis at the expected location of the 
F07 Southeast shelf and displaying the distinct kine- 
matic profile predicted by their orbital model. The sub- 
structure was discovered in the course of an on-going 
Keck/DEIMOS spectroscopic study of the dynamics and 
metallicity of RGB stars in the inner sp heroid and outer 
halo of M31 (see iGilbert et all 120061 . and references 
therein). The portion of the inner spheroid studied 
here appeare d to be relatively undis turbed in earlier star 
count m aps dFerguson et al.l 1200^ ) and radial velocity 
survey s (|Reitzel fc Guhathakurtal 120021 : iKalirai et al.l 
l2006al ). The photometric and spectroscopic data used 
in this analysis are described in §[5J The criteria for se- 
lection of M31 RGB stars are discussed in §[3j The kine- 
matics of the RGB population (first the combined data 
set and then the individual fields) are characterized in §2J 
The spatial trends and general properties of the dynami- 
cally hot spheroid and cold substructure populations are 
discussed in §[5] and §[6j respectively. The physical ori- 
gin of the cold substructure is explored in §0 including 
its likely relation to M31's GSS. The rel evance of the 
newly discovered cold substructure to the iBrown et al.l 
( 2003) discovery of an intermediate-age population in the 
spheroid of M31 is discussed §[H The main conclusions 
of the paper are summarized in §[9] 

2. DATA 

The data set discussed in this paper is drawn from 
photometry and spectroscopy of several fields on/near 
the southeastern minor axis of M31. The locations of the 
fields are shown in Figures Q] and [2j They span a range of 
projected radial distances from M31's center of i? pro j ~ 9 
to 30 kpc (Table [l}. A brief explanation of the data sets 
and reduction is included below. A more detailed dis- 
cussion of the observational strategy and data reduction 
methods employed in our M31 survey can be fo und in 
Guhathakurta et al l (|2006l ) . IKalirai et all (|2006al lbll and 
Gilbert et all (|2006H 



2.1. Photometry 

Photometry and astrometry for the majority of the 
fields analysed in this paper were derived from Mega- 
Cam images in the g' and i' bands obtained with the 
3.6-m Canada-France -Hawaii Telescope ( CFH T) 13 . The 
program SExtractor (jBertin fc Arnoutslll996l ) was used 
for object detection, photometry, and morphological clas- 
sification (via the stellarity parameter). The in- 
strumental g' and i' magnitudes were transformed to 
Johnson-Cousins V and / magnitudes based on obser- 
vations of landolt photometric standard stars (Kalirai 
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Photometry and astrometry for aO, t he outermost field 
discussed in this paper, were derived bv lOstheimerl (2002) 
from images obtained with the Mosaic camera on the 
Kitt Peak National Observatory (KPNO) 14 4-m tele- 
scope in the Washington System M and bands and the 

13 MegaPrime/MegaCam is a joint project of CFHT and 
CEA/DAPNIA, at the Canada- France-Hawaii Telescope which is 
operated by the National Research Council of Canada, the Insti- 
tut National des Science de l'Univers of the Centre National de la 
Recherche Scientifique of France, and the University of Hawaii. 

14 Kitt Peak National Observatory of the National Optical As- 
tronomy Observatory is operated by the Association of Universities 
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Fig. 1. — Sky positions of the fields discussed in this pa- 
per. The blue square represents the position and area of the 
CFHT/MegaCam image (Fig.r2{. The red rectangles approximate 
the size and position angle of the DEIMOS spectroscopic slitmasks. 
The three masks nearest the outer ellipse are in field aO. The re- 
maining m asks are identified in Figure [2] The star count map 
comes fromQBata ct al. ( 20051), and is in standard M31-centric co- 
ordinates 97) . The outer ellipse represents a 55 kpc radius along 
the major axis, with a flattening of 3:5. The major and minor axes 
of M31 are indicated by straight lines. The giant southern stream 
is the obvious overdensity of stars south of M31's center. 

intermediate-width DD051 band. This combination of 
filters allows photometric selection of stars that are likely 
to be M31 red giant branch (RGB) stars rather than 
MW dwarf stars (e.g.. lPalma et al.ll2003l ; lMaiewski et all 
120051 ). The DD051 filter is centered at a wavelength 
of 5150 A with a width of ~ 100 A, and includes the 
surface- gravity sensitive Mg b and MgH stellar absorption 
features which are strong in dwarf stars but weak in RGB 
stars. Based on a star's position in the (M — DD051) 
versus {M—T2) color-color diagram, it is assigned a prob- 
ability of being an M31 RGB star. Johnson-Cousins V 
and I magnitudes were derived from the M and T2 mag- 
ni tudes using the photo metric transformation relations 
in lMaiewski et all (|2000D . Use of DD051 photometry to 
screen for likely M31 RGB stars increases the efficiency 
of the spectroscopic observations, suppressing the num- 
ber of selected MW dwarf stars by a factor of « 3 in aO 
(Guhathakurta et al. 2007, in prep). 

2.2. Spectroscopy 
2.2.1. Slitmask Design and Observations 

Objects in fields covered by the CFHT/MegaCam im- 
ages were selected for Keck/DEIMOS spectroscopy based 
on / magnitude a nd the SExtractor m orphological crite- 
rion stellarity (jKalirai et al.|[2006af ). Objects in field 
aO were selected on the basis of I magnitude and morpho- 
logical criteria (DAOPHOT parameters chi and sharp), 
prioritized according to their probability of being an M31 
RGB star (based on M, T2, and DD051 photometry as 

for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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Fig. 2. — Starcount map derived from CFHT/Mega Cam pho- 
tometry in a single pointing with the 36-CCD mosaic (§ 12.11 . The 
orientation of this map is the same as Figure [l] The size and 
positions of the masks designed from the CFHT/MegaCam pho- 
tometry are shown as red rectangles. The w hite square shows the 
position and approximate orientation of the Brown ct al. (2003) 
deep HST/ACS observations (§[8]l. The three aO masks were based 
on photometry from the KPNO 4-m telescope and are to the south- 
east (bottom left), beyond the limit of this image. There is an 
apparent edge in the density of star counts in the image, running 
from the upper left to the lower righ t an d passing through field 
fl23; this feature will be discussed in § 17.11 



described in § 12. 1[) . Pre-selection of likely M31 RGB stars 
is vital for efficient study of the sparse outer parts of the 
M31 halo. The inner fields (such as the fields drawn 
from the CFHT/MegaCam images) contain a relatively 
high surface density of M31 RGB stars, so the RGB to 
MW dwarf star ratio is high even without DD051-based 
pre-selection of RGB candidates. 

For the purposes of most of the analysis in this paper, 
data are classified according to fields, rather than spec- 
troscopic masks. In general, a "field" refers to the area 
covered by a single CFHT/MegaCam CCD (~ 15'x6f5, 
Fig. [2]); there can be one or more overlapping DEIMOS 
masks (16' x 4') in a single field. For example, masks 
HI 1.1 and Hll_2 are both part of field Hll. There are 
two exceptions to this field/mask scheme: (1) the aO field 
refers to a single Mosaic pointing, which covers a 35' x 35' 
area, and (2) two of the fl30 masks each straddle a couple 
of adjacent MegaCam CCDs, but since they were chosen 
to overlap with the fl30_l mask, they have been labelled 
fl30_2 and fl30_3. 

Fields were observed using the Keck II telescope and 
the DEIMOS instrument with the 1200 line mm -1 grat- 
ing. Most of the fields were observed in Fall 2005 [fl09, 
fll5, fll6, fl23, fl30 (1 mask), and fl35]. The aO masks 
were observed in Fall 2002 and 2004, the Hll masks in 
2004, and the last two of the three fl30 masks in Fall 2006 
(Table [l| . The central wavelength for most masks was 
A7800 A, yielding a spectral coverage of approximately 
AA6450-9150 A. The only exceptions were the a0_l and 
a0_2 masks observed in 2002, which had a central wave- 
length of A8550 A and a spectral coverage of approx- 
imately AA7200-9900 A. The 1200 line mm^ 1 grating 
has a dispersion of 0.33 A pix -1 ; the scale in the spatial 
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direction is 0'.'12 pix -1 , and the effective scale in the dis- 
persion direction is 0'.'21 pix -1 . Slits had a width of 1". 
The spectral resolution for a star observed in typical 0.8" 
seeing conditions is about 1.3 A FWHM. Each mask was 
observed for a total of 1 hour, with the exception of field 
fl09 which was observed for 3 hours. 

2.2.2. Spectroscopic Data Reduction 

The spectra were reduced and analyzed using a mod- 
ified version of the spec2d and spec Id software devel- 
oped by the DEEP2 team at the University of Califor- 
nia, Berkeley 15 ; these routines perform standard spectral 
reduction steps, including flat-fielding, night-sky emis- 
sion line removal, and extraction of the two-dimensional 
spectra. Reduced one-dimensional spectra are cross- 
correlated with a library of template stellar spectra to 
determine the redshift of the object. Each spectrum 
was visually inspected and assigned a quality code based 
on the number and quality of absorption lines. Spectra 
with at least two spectral features (even if one of them 
is marginal) are considered to have secure redshift mea- 
surements. A heliocentric correction is applied to the 
measured radial velocities based on the sky position of 
the mask and the date and time of the observation. The 
heliocentric velocities are not corrected for the changing 
component of solar motion across our fields; our inner- 
most and outermost fields are separated by 1.5° along 
the southeastern minor axis, which corresponds to only 
a 1.6 km s~ x velocity change. 

Spectra in field aO were reduced using the original 
reduction techniques briefly outlined above. Detailed 
discussions of the spectral reduction techniques, qual- 
ity determination, and S/N measurements u sed in our 
survey can be found in iGuhathakurta et al.l (|2006l ) and 
iGilbert et al.1 (|2006l h the typical velocity error in this 
field is 15 km s _1 . Our cross-correlation procedure has 
since been improved and is described below. An off- 
set of +20 km s _1 has been applied to the aO data to 
make them consistent with the results of this new cross- 
correlation procedure. 

Spectra in the remainder of the fields were reduced 
using several improvements to the reduction pipeline. A 
greater number of stellar templates are used for the spec- 
tral cross-correlation, and the template library has been 
expanded to include spectral templates from the Keck 
II telescope's Echellete Spectrograph and Imager (ESI) 
and DEIMOS in addition to the existing Sloan Digital 
Sky Survey (SDSS) spectral templates. The ESI and 
DEIMOS templates were included because they more 
closely match the resolution of the observed spectra. The 
position of the atmospheric A-band in the observed spec- 
trum is used to correct the observed radial velocity for 
imperfect cent ering of the star in the slit (Simon & Geha 
2007, in prep; ISohn et all 12007). The improvement to 
the spectral templates and the A-band correction allows 
us to reduce our velocity measurement errors relative to 
our previous reductions. The median velocity error for 
the data presented in this paper is 4.6 km s _1 , estimated 
from the cross-correlation output routine and confirmed 
by repeat measurements of individual stars on overlap- 
ping DEIMOS masks. 

15 http : //astron . berkeley . edu/~cooper/deep/spec2d/primer . 
http : //astron .berkeley . edu/~cooper/deep/specld/primer .html 



Spectroscopic data in fiel ds Hll and aO have b een pre- 
sente d in previous papers dKalirai et al.ll2006al: Brown 
et al. 120031 . 12006af ). The radial velocity sample for Hll 
presented in this paper contains « 50% more M31 RGB 
stars than the previously published sample, due to the 
recent recovery of spectra from two CCDs that were not 
included in these earlier papers and the improvements in 
the data reduction process described above. 

3. SELECTING A SAMPLE OF M31 RED GIANTS 

The largest source of contaminants in our spectroscopic 
survey are foreground MW dwarf stars. Background 
galaxies are easy to identify and remove from the sam- 
ple on the basis of their spectra and redshifts. However, 
the radial velocity distribution of MW dwarf stars over- 
laps that of M31 RGB stars, making identifying individ- 
ual stars as M31 red giants or MW dwarf s problematic. 
We us e the diagnostic method detailed in IGilbert et al.l 
(120061 ) to separate M31 RGB stars from MW dwarf star 
contaminants. The method uses empirical probability 
distribution functions to estimate the likelihood a given 
star is an M31 red giant based on five photometric and 
spectroscopic diagnostics: 

• The radial velocity of the star. 

• Photometry in the M, T%, and (surface-gravity sen- 
sitive) DD051 bands 

• The measured equivalent width of the Na I doublet 
at 8190A combined with the (V — I)o color of the 
star. 

• The position of the star in an (I, V — I) color mag- 
nitude diagram with respect to theoretical RGB 
isochrones. 

• A comparison of the star's photometric vs. spec- 
troscopic metallicity estimates. 

The DD051 diagnostic is only used for field aO, which 
is the only field in the present work for which DD051 
photometry is available (§ 12. ip . The likelihoods for each 
diagnostic are combined in a weighted average for each 
star to determine the overall likelihood, (Li), the star is 
an M31 RGB or MW dwarf star (§ A.l). Based on the 
overall likelihood, each star is identified as either a se- 
cure M31 RGB star ((L t )> 0.5, or > 3x more likely 
to be an M31 RGB star than an MW dwarf) or se- 
cure MW dwarf star ((Li)< —0.5), or a marginal M31 
RGB star (0 <(L t )< 0.5 ) or marginal MW dwarf star 
(-0.5 <(Li)< 0). 

An advantage of the diagnostic method for the present 
analysis is the ability to select a sample that is chosen 
independently of radial velocity. Since radial velocity is 
only one of a number of diagnostics which are used to de- 
termine the nature of an individual star, the likelihood 
method (even with the inclusion of the radial velocity di- 
agnostic) presents a significant improvement over the use 
of velocity cuts to select samples for kinematical analy- 
sis by reducing the sensitivity of the sample to velocity. 
However, by using the likelihood method without the ra- 
dial velocity diagnostic (resulting in overall likelihoods 
(Li)S), we are able to select a sample of M31 red giants 
html that is completely independent of the radial velocities of 
the stars. 
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TABLE 1 

Details of Spectroscopic Observations and Basic Results. 



Mask 


I-' vr\ i ac'i- &c\ 


Pointing center: 


PA 




=ft Sri 


44- n f TVTQI 




Radius 


.12000 


^J2000 


[ hi OI IN ) 


Ubs. (U1J 


targets a 


otars ' 




(kpc) 


\ ' " / 


\ ■ ■ / 










fl09 1 


9 


00:45:46.75 


+40:56:53.8 


23.90 


900^ Alio- 9Q 


208 


169 


mill 


12 


00:46:21.02 


+40:41:31.3 


21.0 


2004 Sep 20 


139 


89 


Hll_2 


12 


00:46:21.02 


+40:41:31.3 


—21.0 


2004 Sep 20 


138 


88 


f 116.1 


13 


00:46:54.53 


+40:41:29.5 


22.60 


2005 Aug 28 


199 


119 


f 115.1 


15 


00:47:32.71 


+40:42:00.9 


-20.0 


2005 Aug 28 


191 


111 


fl23_l 


18 


00:48:05.57 


+40:27:16.3 


-20.0 


2005 Aug 28 


171 


101 


fl35.1 


18 


00:46:24.88 


+40:11:35.5 


-27.0 


2005 Aug 29 


146 


99 


fl30.1 


22 


00:49:11.97 


+40:11:45.3 


-20.0 


2005 Aug 28 


108 


52 


fl30.2 


23 


00:49:37.49 


+40:16:07.0 


90.0 


2006 Nov 21 


115 


13 


fl30.3 


20 


00:48:34.59 


+40:16:07.0 


90.0 


2006 Nov 22 


124 


11 


aO.l 


31 


00:51:51.32 


+39:50:21.4 


-17.9 


2002 Aug 16 


89 


25 


a0.2 


31 


00:51:29.59 


+39:44:00.8 


90.0 


2002 Oct 12 


89 


32 


a0_3 


29 


00:51:50.46 


+40:07:00.9 


0.0 


2004 Jun 17 


90 


26 



a A number of targets were observed on two different masks. Therefore, the total number of unique 
science targets/M31 RGB stars in fields Hll, fl30 and aO is less than the reported number. There 
are 18 M31 RGB stars with duplicate observations (2 in Hll, 8 in fl30, and 8 in aO), thus the 
total number of unique M31 RGB stars is 1013. b The number of M31 RGB stars is defined as the 
numb er of stars that are identified as secure and marginal M31 RGB stars by the Gilbert et al. 
(2006) diagnostic method, wihout the use of the radial velocity diagnostic ({-£<i)y> 0, §[3}. 



Fig. [3] presents the radial velocity distribution of 
stars selected as M31 red giants based on the diagnos- 
tic method with (shaded histograms) and without (thick 
open histograms) radial velocity included, for multiple 
(Li) and (£i)y thresholds. For reference, the radial ve- 
locity distribution of all stars with successful velocity 
measurements is also shown (thin open histograms); the 
MW dwarf star contaminants form the secondary peak 
centered at Vh c \ « —50 km s _1 . The M31 RGB distribu- 
tions are similar, with the sample that includes the ra- 
dial velocity diagnostic showing a systematic deficiency 
of stars at radial velocities near km s _1 . The effect 
of the radial velocity diagnostic on the overall likelihood 
((Li) and (£i)y) distributions is discussed in § A.l. 

The M31 RGB samples identified by their (Li) values 
have a minimal amount of MW dwarf star contamina- 
tion, but are also kinematically biased against stars with 
velocities near km s" 1 (§A.l). The M31 RGB sam- 
ples identified by their (Li)y values have a slightly larger 
amount of dwarf contamination (particularly evident in 
the bottom panel of Figure [3]), but the underlying M31 
RGB population is kinematically unbiased. 

The RGB sample used in this paper is defined as stars 
that are designated as secure and marginal M31 red gi- 
ants by the diagnostic method, with the radial velocity 
likelihood not included in the calculation of a star's over- 
all likelihood of being an M31 RGB star (i.e., (£,)y> 0). 
The number of M31 RGB stars in each field is listed in 
Table [T] The (Li)y> threshold maximizes the com- 
pleteness of the underlying, kinematically unbiased M31 
RGB population, but introduces an overall MW dwarf 
star contamination of 5% to the sample (§ A. 2). The con- 
tamination is largely constrained to Vhd > —150 km s _1 
due to the velocity distribution of MW dwarf stars, and 
its effect on the measured parameters of the M31 RGB 
sample is quantified in § A.3. 

4. STELLAR KINEMATICS IN M31'S SOUTHEAST 
MINOR-AXIS FIELDS 
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Fig. 3. — The radial velocity distributions of samples with 
(shaded/ dotted histograms) and without (thick open histograms) 
the radial velocity diagnostic included in the overall likelihood cal- 
culation. The radial velocity distribution of all stars with suc- 
cessful radial velocity measurements is shown for comparison (thin 
open histograms); the MW dwarf star contaminants form the 
peak at % e ] fa —50 km s —1 . Top: Stars designated as secure 
([(Li),(Li)y]> 0.5) M31 red giants only. Middle: Stars designated 
as marginal and secure M31 red giants ([(Li),{Li)u]> 0). Bot- 
tom: Stars designated as marginal MW dwarfs, marginal M31 red 
giants, and secure M31 red giants ([{L;),{I/;)y]> —0.5). The dis- 
tributions are similar, but there is a deficiency of stars at velocities 
near km s _1 in the sample that includes radial velocity as a 
diagnostic. 

The data presented in this paper span a range in pro- 
jected radial distance from the center of M31 of 9 to 
30 kpc, along the southeastern minor axis. We refer 
to this region as the "inner spheroid," even though it 
has traditionally been referred to as the "halo" of M31. 
This region departs from the classical picture of a stel- 
lar halo formed from observations of the Milky Way: 
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M31's inner sphe roid is about 10 x more metal-rich than 
the M W's halo ([Mould fc Kristianl U986I : iDurrell et~aTl 
l2004f l a nd follows a de Vaucouleurs r 1 / 4 s urface density 
profil e (|Pritchet fc van den Berghl 119941 ; IDurrell et al.l 
120041 ) . while the surface density profile of the MW halo 
follows an r~ 2 power law. In other words, M31's in- 
ner spheroid appears to be a continuation of its central 
bulge. Recent observations have discovered an outer stel- 
lar "halo" in M31 which is relatively metal-poor (Kalirai 
et al. [2006b; C hapman et al.ll2006l ). has a surface density 
profi le that follow s an ~ r~ 2 power law f Guhathakurta 
et al. Ill I: llrwin et al. 2005) , and has b een detected out 
to i? proj = 165 kpc (|Gilbert et al.|[2006h . 

These observations imply that the spheroid of M31 
has two components: an inner, de Vaucoulcurs profile 
spheroid, and an outer, power-law profile halo. We thus 
use the term inner spheroid to distinguish the region 
-Rproj ~ 9 - 30 kpc from the canonical central bulge and 
the newly discovered stellar halo of M31. The outer 
limit of this region is well-defined; a break in the sur- 
face brightness profile of M31 RGB stars has been ob- 
served at -Rnrni ~ 20 — 30 kpc ([Guhathakurta et"aLll2005l : 
llrwin et al.[|2005l ), and there is observational evidence 
that the crossover between the predominantly metal- 
rich population of the inner spheroid and the predom- 
inantly metal-poor population of the outer halo occurs 
at i? proj ~ 30 kpc ([Kalirai et alj|2006bl ). The inner limit 
of this region is arbitrary, as the relationship between 
this component and the central bulge of M31 is not yet 
clear. 

The rest of this section characterizes the line of sight 
velocity distribution of stars in the inner spheroid of 
M31 through maximum-likelihood fits of Gaussians to 
the combined data set f£j |4.1[) and to individual fields 
(§ 14. 2|) . Gaussians provide a convenient means of fit- 
ting for multiple kinematical components in the data 
and characterizing their mean velocity and velocity dis- 
persion. The true shape of the velocity distribution of 
a structural component in M31 is likely to be different 
from a pure Gaussian. However, given the limited sam- 
ple size and the absence of any specific physical model, 
the choice of Gaussians seems appropriate. 

4.1. Maximum- Likelihood Fits to the Velocity 
Distribution of the Combined Data Set 

Figure [4] shows the combined radial velocity distribu- 
tion of M31 RGB stars from all eight fields along the 
minor axis, ranging from 9 to 30 kpc in projected radial 
distance from the center of M31. Fits to the radial veloc- 
ity distribution were made using a maximum-likelihood 
technique; the best-fit single (a) and double (6) Gaus- 
sians are displayed in Figure [U A reduced x 2 analysis 
rules out the single-Gaussian fit, as the probability is 
that the observed radial velocities were drawn from 
the best-fit distribution. The observed velocity distribu- 
tion is well fit by a sum of two Gaussians (solid curve, 
panel b of Fig. 0|, composed of a wide Gaussian (dashed 
curve) centered at (v) sph = -287.2±*;° km s -1 , with a 
width of cr^ ph = 128.9lg'g km s _1 , and a narrow Gaussian 
(dotted curve) centered at (v) suh = -285 AtHi km s" 1 
with a width of <7™ b =42.2^43 km s _1 , which comprises 
19^8% of the total population (quoted errors represent 
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Fig. 4.— The radial velocity distribution of the M31 RGB 
inner spheroid population. A maximum-likelihood analysis was 
used to fit an analytic function to the distribution. (a) The 
best fit single Gaussian has parameters (ri) s P h = —287 km s — 1 and 

<Tv ph =117 km s — 1 . A Xr e d * es * ru l es ou t the single-Gaussian fit 
at a very high confidence level. (ft) The best constrained double- 
Gaussian fit (Table [2} is shown as a solid curve, with the narrow 
component and wide components displayed separately as dotted 
and dashed curves, respectively. 

the 90% confidence limits from the maximum-likelihood 
analysis). Due to the MW dwarf star contaminants in the 
M31 RGB sample (§EJ), the true (v) sph value of the wide 
Gaussian component is 15 to 20 km s" 1 more negative 
than the best-fit value (§ A. 3), making it consistent with 
the systemic velocity of M31 (v sys — —300 km s _1 ). The 
kinematically hot component corresponds to the inner 
spheroid of M31 (quantities related to this component are 
denoted with the subscript "sph"), while the kinemati- 
cally cold component corresponds to substructure in the 
inner spheroid (denoted with the subscript "sub"); the 
discussion of these components and their properties will 
be deferred to §[5] and §[51 respectively. The wider of the 
two Gaussian components in the double Gaussian fit to 
the combined data set is hereafter referred to as G sph (v). 

Figure [5] shows error estimates from the maximum- 
likelihood analysis (in the form of A\ 2 = X 2 ~ Xmin 
curves) for the five double-Gaussian parameters. The 
best-fit value of each parameter is marked as well as 
the 90% confidence limits. The A% 2 curves have a deep 
minimum for all five parameters, an indication that the 
double-Gaussian model is a good description of the ob- 
served radial velocity distribution of inner spheroid stars. 

4.2. Maximum-Likelihood Fits to the Velocity 
Distributions of Individual Fields 

As discussed in the previous section, the combined 
data set shows definite evidence of a kinematically cold 
component centered at about —300 km s _1 that com- 
prises a significant fraction (19%) of the total popula- 
tion. We next investigate which of the fields in our data 
set are the main contributors to this cold population. 
Figure [6] shows velocity histograms for each of the eight 
fields analyzed in this paper. While most of the fields 
display hints of substructure — in the form of one or 
more small possible peaks in their velocity distribution 
that may be marginally significant relative to the (sub- 
stantial) Poisson noise — we are specifically interested in 
judging each field's contribution to the cold component 
at ~ —300 km s _1 . For this purpose, we compare the 
data in each field to three sets of Gaussian fits. We de- 
scribe each of the fits below, summarize the results of the 
fits, and then list the quantitative details of the fits for 
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TABLE 2 

Radial Velocity Distributions: Best Fit Gaussian Parameters. 



Field 



Best fit Gaussian Parameters' 1 
Cold Component Hot Spheroid b 



(v) 



s u hi 



{v) a P h 



sph 
rr * 



Fraction 

N suh /N to t 



All fields 


-285.4^.8 


fl09 




Hll 


-294.3li^ 


fll6 


-309.4j;^;g 


fll5 




fl23 


-279.4t^i 


fl35 




fl30 




aO 





42.2; 
55.5; 

51.2 
c lO.i 



-12.5 
-14.3 



15.6 
12.7 
+24.4 
15.0 



H6.9 
-5.0 



-287.2 
-274.5 



+8.0 

-7.7 
15.4 
15.3 

-287.2 

-287.2 

-270.9±i|;« 

-287.2 

oro.i_ 21 3 
s +19.3 
-19.2 
,+25.2 
25.2 



-259. 8Z 

-299. 2 H 



128.9 
120.7 



+7.7 

-6.9 
11.7 
10.1 

128.9 
128.9 
120 1+ 14 - 4 
128.9 

131.5±i|"l 

10, c+29:5 

101.0_ 21 g 



0.19 



-0.09 
-0.08 



+0.16 
-0.16 

-1-0.22 



0.31 



+0.11 



a A double Gaussian fi t is presented for the combined sample, as it is a poor fit 
to a single Gaussian (§ 14. 1^ . Constrained double Gaussian fits are presented for 
three of the fields (Hll, fll6, and fl23), with the wider component held fixed 
(adopting the fit to the combined sample). Single Ga ussia n fits are presented 
for the remaining five fields. The reader is referred to § 14. 21 for details of the fits 
to individual fields. Errors quoted represent the 90% confidence limits from the 
maximum-likelihood analysis. 15 The M31 RGB sample used in this analysis was 
chosen to ensure a high degree of completeness, and thus suffers from some MW 
dwarf contamination (§A.2). The MW dwarf star contaminants are largely at 
v hcl > —150 km s —1 , and cause the best-fit (i)} sph values to be biased towards 
more positive velocities. The true (i;) sph values of the M31 RGB population 

are 15 to 20 km s — 1 more negative than listed here. The effect on ov P is 
negligible (§ A.3). c The median velocity error for the stars in fl23 is 4.6 km s _1 . 
The estimated intrinsic velocity dispersion of the cold component in fl23 after 
accounting for velocity measurement error is 9.5 km s — . 
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Fig. 5. — Results of the maximum-likelihood analysis for the 
double-Gaussian fit to the combined M31 RGB inner spheroid sam- 
ple. The optimal value of each parameter is marked with an arrow, 
and the 90% confidence limits from the maximum-likelihood analy- 
sis are marked with dashed lines. The upper limit of the y-axis rep- 
resents the 99% confidence limits. The parameter A.\ 2 = X 2- Xmin 
is plotted as a function of (a) mean velocity of the cold substruc- 
ture component (narrow Gaussian), (u) sub , (f>) velocity dispersion 
of the cold component, crj ub , (c) mean velocity of the M31 in- 
ner spheroid (wide Gaussian), (?;} sph , (d) velocity dispersion of 
the M31 inner spheroid, av ph ! and (e) fraction of the total M31 
RGB population in the cold component, N au ^/Ntot- The horizon- 
tal arrow in panel (c) represents the correction to the (u) aph value 
necessary to offset the bias caused by MW dwarf contamination at 
v hcl > -150 km s" 1 (§A.3). 
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Fig. 6. — Velocity histograms for each of the individual fields 
with best-fit Gaussians overlaid. Fields that did not sh ow c lear 
evidence of substructure (fl09, fll5, fl30, fl35 and aO, § 14.21 1 are 
shown with the Gaussian component from the double-Gaussian fit 
to the full sample (G Bph (v), solid curves), as well as the best- 
fit single Gaussian to the individual field (dotted curves). Fields 
with evidence of substructure (Hll, fll6, fl23), are shown with 
both their best-fit single (dotted curves) and double (solid curves) 
Gaussians. For the constrained double-Gaussian fits, both the nar- 
row and broad (G Bph (v)) Gaussian components are shown (dashed 
curves) scaled to their relative contributions. 



each field. 

First, the radial velocity distribution in each field is 
compared to the Gaussian G sph (v) defined by (v) sph = 
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-287.2 km s~\ <r s v ph = 128.9 km s" 1 (§ OJ using 
the reduced x 2 statistic (x 2 od ). Fields f!09, fl 15, and 
aO are consistent with being drawn from G sph (v), and 
so are ruled out as significant contributors to the ~ 
—300 km s _1 cold component. The rest of the fields are 
at least marginally inconsistent with being drawn from 
G sph (v). 

Second, a maximum-likelihood single-Gaussian fit is 
performed on the radial velocity distribution in each field 
and compared to the data. Fields fl23 and fl35 are in- 
consistent with their respective best-fit single Gaussians 
(based on the x 2 ed statistic) and are therefore suspected 
to contain substructure. The remaining fields (fl09, Hll, 
fl 16, fl 15, fl30, and aO) are consistent with their best- 
fit single Gaussians. For fields fl09, fll5, fl30, and aO, 
the best- fit Gaussians are consistent with G aph (v), and 
they are thus ruled out as significant contributors to the 
~ —300 km s _1 cold component. Fields fll6 and Hll are 
suspected to contain substructure because their best-fit 
Gaussians are significantly narrower than G sph (v). Pre- 
vious kinematic studies of M31's inner spheroid, includ- 
ing this one, have found a ubiquitously broad distribu- 
tion of radial velocities (v\ lc \ ~ to w —600 km s _1 ; 
§ [5]). Thus, the anomalously narrow single Gaussian fits 
in fields fll6 and Hll cause us to suspect them of being 
contributors to the ~ —300 km s^ 1 cold component. 

Third, we carry out a maximum-likelihood fit to all 
fields using a constrained double Gaussian, with G sph (i>) 
as the fixed wide Gaussian component. Fields Hll, fll6, 
and fl23 are well fit by a constrained double Gaussian 
(based on the x 2 ed statistic and well-defined minima for 
the variable Gaussian parameters). These three fields 
are significant contributors to the ~ —300 km s _1 cold 
component, and the constrained double-Gaussian fit is 
adopted as the preferred fit in the subsequent discussion. 
The best-fit cold component in the constrained double- 
Gaussian fit to fl30 is centered at ~ —50 km s _1 , but this 
likely represents residual contamination by MW dwarf 
stars (§EJ §A). The remaining fields (fl09, fll5, fl35, 
and aO) are poor fits to a constrained double-Gaussian 
in that the Gaussian parameters do not have well-defined 
minima. 

In summary: 

• Three fields, Hll, fl 16, and fl23 are identified as 
significant contributors to the ~ —300 km s _1 cold 
component. 

• Although field fl35 shows evidence of substructure, 
it is a poor fit to a constrained double- Gaussian and 
is not a definite contributor to the ~ —300 km s _1 
cold component. (However, the fit may be con- 
fused by the presence of multiple cold components; 

mm) 

• Four fields, fl09, fll5, fl30, and aO, are not sig- 
nificant contributors to the ~ —300 km s" 1 cold 
component. 

Each field is discussed individually below, and Table [2] 
summarizes the preferred fits (single or double-Gaussian) 
to the velocity distributions in each field. 
Field fl09: The data in this field are consistent with be- 
ing drawn from the Gaussian G aph (v) f§ l4.ip . The best-fit 
single Gaussian to the data in this field has parameters 



(v)= -274.5 km s _1 and a = 120.7 km s _1 , and the 
data are also consistent with being drawn from this dis- 
tribution. 

Field Hll: The probability the data in this field are 
drawn from the Gaussian G sph (v) is P < 1%, thus the 
data are inconsistent with this distribution. The best-fit 
single Gaussian to the data in this field has parameters 
(v)= -291.1 ± 11.6 km s" 1 and a = 106.2j^ km s _1 , 
and the data are consistent with being drawn from this 
distribution. However, the best-fit value of a in this 
field is anomalously low compared to the value of tr^ ph 
determined from the double-Gaussian fit to the com- 
bined data set: the two values are inconsistent at the 
~ 3.5cr level. This suggests that there is a kinemat- 
ically cold population in this field, and a comparison 
of the data to the constrained double-Gaussian fit to 
this field ((w) sub = -294.3 km s"\ a s v uh = 55.5 km s _1 , 
As U b/-/Vtot = 0.44) returns a x^ed ^ na * ^ s significantly 
smaller than the x 2 ed of the single Gaussian fit. 
Field fll6: The data in this field are inconsistent with 
being drawn from the Gaussian G sph (v). The best- fit 
single Gaussian to the data in this field has parameters 
(v)= -292.9±{{i km s" 1 and a = 97±rl km s"\ and 
the data are consistent with being drawn from this dis- 
tribution. As in field Hll, the best-fit value of a is in- 
consistent with cr^ ph from the combined data set, at the 
~ 4.2a level. A comparison of the data to the constrained 
double-Gaussian fit to this field ((u) sub = -309.4 km s~\ 
a s v uh = 51.2 km s~\ N suh /N tot = 0.44) returns a signifi- 
cantly smaller x 2 cd than that of the single-Gaussian fit. 
Field fll5: The data in this field are consistent with 
being drawn from the Gaussian G sp (v), as well as the 
best-fit single Gaussian ((u) sph = -270.9 km s~\ < ph = 
120.1 km s" 1 ). 

Field f 123: The probability that the data in this field are 
drawn from the Gaussian G sph (v) is P « 1%. The prob- 
ability that the data are drawn from the best-fit single 
Gaussian ((w) sph = -270.9 km s _1 , cr sph = 120.1 km s" 1 ) 
in this field is also « 1%. The data in this field 
are strongly inconsistent with being drawn from any 
single Gaussian, but they are consistent with the con- 
strained double-Gaussian fit ((u) sub = —279.4 km s _1 , 
of ub = 10.6 km a -1 , N suh /N tot = 0.31). The median 
velocity error for the stars in field fl23 is 4.6 km s _1 , 
therefore the intrinsic velocity dispersion of the cold com- 
ponent in this field is estimated to be 9.5 km s _1 . 
Field fl35: The data in this field are inconsistent with 
being drawn from the Gaussian G sph (v) (P ~ 1%). 
The best-fit single Gaussian to the data in this field 
has parameters (v) sph = —315.1 km s _1 and <Tv ph — 
127.8 km s _1 , but the data are also inconsistent with 
being drawn from this distribution (P < 1%). How- 
ever, the maximum-likelihood analysis was unable to 
constrain a double-Gaussian fit to any reasonable degree 
of certainty — i.e., the error estimates on the parameters 
show no strong global minima. Field fl35 is therefore 
treated as a field without a definite detection of substruc- 
ture, although in § !7.1.3l we discuss the possible presence 
of multiple kinematically-cold components in this field. 
Field flSO: The data in this field are inconsistent with 
being drawn from the Gaussian G sph (v) (P < 1%), but 
are consistent with being drawn from the best-fit single 
Gaussian ((v) Bph = -259.8 kms" 1 , a sph = 131.5 kms" 1 ). 
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The best-fit single Gaussian for this field has parameters 
which are consistent with G sph (v) (Table [2]). 
Field aO: The data in this field are consistent with be- 
ing drawn from the Gaussian G sph (v) as well as the 
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best-fit single Gaussian ((w) sph = 
131.5 km s- 1 ). 
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We now have a characterization of the kinematical 
properties of the combined data set and the individual 
fields, and have identified the fields which significantly 
contribute to the ~ —300 km s _1 cold component dis- 
covered in the combined data set (Hll, fll6, and fl23). 
The number of stars in fields Hll, fl 16, and fl23 which 
are associated with the cold component comprise 17% of 
the total number of M31 RGB stars in the sample, con- 
firming that these fields are the primary contributors to 
the cold component. Next, we discuss the trends in the 
properties of the dynamically hot and cold populations 
(§[5] and §|6l respectively), followed by possible physical 
interpretations of the cold component (§[7J). 

5. VELOCITY DISPERSION OF M31'S VIRIALIZED 
INNER SPHEROID 

As discussed in §[H the radial velocity distribution of 
the combined data set contains a significant cold compo- 
nent at Uhoi ~ —300 km s _1 . An analysis of the kinemat- 
ical profile of the individual fields shows that this is due 
to the presence of a significant amount of substructure 
in three of the fields (§ 14. 2p . The kinematically hot com- 
ponent in the double Gaussian fit to the combined data 
set presumably represents the underlying virialized inner 
spheroid of M31. The velocity dispersion of the spheroid 
based on a maximum-likelihood double Gaussian fit to 
the combined data set (1013 M31 RGB stars) is < ph = 
128.9 km s -1 . 

We have combined all of the fields to get a more 
robust estimate of the velocity dispersion of the inner 
spheroid, since the individual fields suffer from small 
number statistics. However, the dynamical quantities 
((i>) sph , <7y ph ) may have radial dependencies. The veloc- 
ity dispersion of sphero ids is expected to decrease with 
increasing ra dius (e.g.. iNavarro "et~aT1fT99l. l200l Ma- 
mon & Lokas 2005t iDekel et al J 120051 ). Rotationof the 
spheroid or the tangential motion of M31 could result 
in a radial dependency of (i>) sph , but since our fields are 
mostly aligned along the minor axis and cover a relatively 
small radial range, we do not expect to be sensitive to 
these effects. 

To test for a dependency of the dynamical quantities 
on radius, we analyze the single Gaussian fits to the ve- 
locity distributions in each of the five fields which do 
not show clear evidence of substructure as described in 
§112] (Table (2). The Ax 2 error estimates for the best- 
fit single Gaussian parameters in these fields are shown 
in Figure [7J The best-fit (v) sph values in each field 
are largely consistent with the systemic velocity of M31 
(v sys — —300 km s -1 ) once the effect of MW dwarf 
star contamination is taken into account, with the ex- 
ception of fl35, whose (v) sph is significantly more nega- 
tive. Our data show no evidence for a decreasing a^ ph 
with increasing radius. However, the fields discussed in 
this paper only span a range in projected radial dis- 
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Fig. 7. — Results of the maximum-likelihood analysis for the 
single-Gaussian fits to individual fields that do not show clear ev- 
idence of the ~ —300 km s _1 cold component in their radial ve- 
locity distributions fij |4.2t . The panels show Ax 2 for the mean 
velocity (-u) sph (left) and velocity dispersion ov ph (Tight) of M31 
RGB stars in fields (a) fl09, (b) fll5, (c) fl35, (d) fl30, and (e) 

aO. As in Figure \E\ the optimal values of (?;) sph and ov P are 
marked by arrows, and the 90% confidence limits are shown as 
dashed lines. The horizontal arrow in each of the left panels repre- 
sents the correction (magnitude and direction) to the (i>) s P h value 
necessary to offset the bias caused by MW dwarf contamination at 
v hcl > -150 km s" 1 (§ A. 3). 



to the size of the total M31 spheroid (inner spheroid and 
outer halo), which h a s been shown to extend to 165 kpc 
(Gil bert et aT][20M ). iBattaglia et~aH <|2005) use a sam- 
ple of 240 Galactic halo objects to determine the radial 
velocity dispersion of the Milky Way halo, and find that 
it has an almost constant value of 120 km s _1 out to 
30 kpc, beyond which it decreases with increasing radial 
distance, declining to 50 km s _1 at 120 kpc. 

There have been a few previous measurements of 
the velocity dispersion based on M31 spheroid stars. 
The closest a nalog to the present study is Reitzel & 
Guhathakurta fl2002j), who fit 80 candidate M31 RGB 
stars at i? pro j = 19 kpc on the southeastern minor 
axis with a combination of the Galactic standard model 
and a wide Gaussian. They found cr^ ph ^ 150 km s _1 
for the M31 velocity dispersion, with the number of 
M31 RGB stars estimated to be 43% of the population. 
iGuhathakurta et~al . (2006) studied a field on the giant 



southern stream at i? pr oj = 33 kpc and found a veloc- 
ity dispersion of tr^ ph = 65^1 km s_1 f° r the underlying 
spheroid based on a sample of « 21 stars. However, if 3 
likely RGB stars with % e i > —150 km s _1 are included, 
the estimated velocity dispersion in this field increases to 



sph_ 



116t^ km s~\ 



tance of i? D 



9 to 30 kpc. This is small compared 



Chapman et al. (|2006f ) have measured a mean velocity 
dispersion of cr^ ph = 126 km s^ 1 for the inner spheroid 
of M31 using ~ 800 RGB stars in multiple fields sur- 
rounding M31. They determine that the spheroid has a 
central velocity dispersion of 152 km s" 1 which decreases 
by —0.9 km s _1 kpc -1 out to i? pro j ~ 70 kpc. Many of 
their fields are near M31's major axis and have signifi- 
cant contamination from the extended rotating compo- 
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nent identified as the extended disk in llbata et al.1 (2005). 
They isolate a sample of M31 spheroid stars by removing 
all stars within 160 km s _1 of the disk velocity in each 
field, and by removing all stars with Uhci > —160 km s _1 
(MW dwarf star contaminants). This "windowing" tech- 
nique leaves them with a sample of M31 spheroid stars 
that is significantly incomplete, but is largely uncontami- 
nated by M31's extended disk o r MW dwarf stars. Based 
on the iChapman et al.1 (|2006l ) result, we would expect 
to measure a velocity dispersion of 146 km s _1 in our 
innermost field, fl09, which would decrease to a veloc- 
ity dispersion of 125 km s _1 in our outermost field, aO. 
The predicted velocity dispersion of 146 km s _1 in field 
fl09 exceeds the 90% confidence limits of the maximum- 
likelihood fit and is just within the 99% confidence limit, 
and the data presented in this paper show no evidence 
of a strong trend in cr^ ph with radius. 

The previous measurements of the stellar velocity dis- 
persion of the M31 spheroid relied either on samples that 
were chosen on the basis of radial velocity cuts or on a 
statistical fit to the combined M31 RGB and MW dwarf 
populations. Our measurement of the velocity disper- 
sion of M31's spheroid is unique in that it is based on 
a sample of spectroscopically confirmed M31 RGB stars 
that were chosen without the use of radial velocity (§[3]). 
Our method also allows us to quantify and correct for 
the effect of MW dwarf star contamination (§A.2 and 
§A.3). 

6. PROPERTIES OF THE MINOR-AXIS 
SUBSTRUCTURE 

6.1. Spatial Trends: Kinematics and Structure 

As discussed in §|4~2j fields Hll, fll6, and fl23 show 
evidence of the cold component at ~ —300 km s _1 in 
their radial velocity distributions (Fig. [6]) , and are well fit 
by a sum of two Gaussians with G sph (v) as the fixed wide 
Gaussian component (Table [2]). The A% 2 error estimates 
for the free parameters ((v) suh , cr™ b , and N su b/N to t) are 
shown in Figure [8] for each field. 

If the three fields are considered together, a pattern 
emerges. Both the velocity dispersion and the fraction 
of stars in the cold component decrease with increas- 
ing radial distance. The cold component in fields Hll 
and fl 16 (at i? pro j = 12 and 13 kpc, respectively) is sig- 
nificantly wider than the cold component in field fl23 
(i? pro j = 18 kpc). The cold component also appears to 
be more dominant by number (surface brightness) over 
the hot component in fields Hll and fll6 than in field 
fl23. 

Figure [9] presents the velocities of the M31 RGB stars 
as a function of their distance along the major and mi- 
nor axes of M31. The cold component can be seen as a 
triangular-shaped feature that narrows to a sharp point 
as the distance along the minor axis increases. The fields 
that overlap the triangular-shaped feature in minor axis 
distance are (in order of increasing distance along the mi- 
nor axis) Hll, fll6, fll5, fl35, and fl23. Stars in fields 
Hll, fll6 and fl23 that are within the area denoted by 
the dotted line are shown as red crosses in the bottom 
panel. In the fields in which it is observed, the cold com- 
ponent is spread evenly along the direction of the major 
axis, and is centered at «hoi^ —285 km s . The major- 
ity of the fields overlap in position along the major axis; 
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Fig. 8. — Results of the maximum-likelihood analysis for the 
narrow Gaussian parameters from the constrained double-Gaussian 
fits to fields (a) Hll, (6) fll6, and (c) fl23. The panels show 
Ax 2 as a function of (left) mean velocity (i>) sub , (middle) velocity 
dispersion cr Buh , and (right) fraction of stars in the cold component, 
Aaub/^tot- The parameters of the wide Gaussian component in 
these fits have been fixed at the values of the Gaussian G Bph (v) 
(§ 14,11 1, As in Figure [5] the optimal values of each parameter are 
marked by arrows, and the 90% confidence limits are shown as 
dashed lines. The velocity dispersion decreases with increasing 
radial distance from the center of M31. 

the exception is field fl35, which is the isolated set of 
points at large major axis distances (—0.51° to —0.36°, 
or —7.0 kpc to —4.9 kpc) in the bottom panel of Figure[H 
Field fl23 extends from -1.18° to -1.39° (16.1 kpc to 
19.0 kpc) along the minor axis; the tip of the feature is 
seen in this field. This field also has the coldest substruc- 
ture in the radial velocity histograms (Fig. [51 Table [5]) . 
Field Hll brackets the other edge of the feature along 
the minor axis. Field Hll has a range in minor axis dis- 
tance of -0.78° to -0.99° (10.7 kpc to 13.5 kpc), and 
has a large degree of overlap with field fl 16 in both their 
minor and major axis distance ranges. The fits to the 
cold component in these two fields return similar <7^ ub 
and N snb /N tot estimates (Table [2]). 

Of the fields in which substructure was not clearly de- 
tected kinematically (§ I4.2[) . fields fl35 and fll5 are both 
within the minor axis distance range spanned by the ob- 
served substructure (dotted triangular region in upper 
panel of Fig. [5]). A two-sided Kolmogorov-Smirnov (KS) 
test finds that the radial velocity distribution of field fll5 
is consistent with the radial velocity distribution of its 
closest neighbor, field fll6 (Fig. [5]), although it is in- 
consistent with the best constrained double Gaussian fit 
to field fll6. Although the radial velocity distribution 
of field fl35 is not well-fit by a double Gaussian (§ 14. 2|) . 
there is a concentration of stars near Vhci ~ —300 km s _1 
in its radial velocity distribution (Fig. [6]). We will discuss 
these fields further in the context of the physical inter- 
pretation of the cold substructure (§ I7.1|) . Fields fl30 and 
aO are at larger minor axis distances than the tip of the 
feature, and field fl09 is interior to the feature. The fact 
that substructure is not detected in these fields is con- 
sistent with our favored physical interpretation of the 
—300 km s~ x cold component, which will be discussed 
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Fig. 9. — Distribution of M31 RGB stars in velocity vs. distance 
along the minor (top) and major (bottom) axes of M31. The range 
of minor (major) axis distances of stars in each field are shown in 
the top (bottom) panel. Top: The cold component is visible as a 
triangular feature that starts at —0.8° (10.9 kpc) and narrows to 
a point at —1.35° (18.4 kpc) along the minor axis. This feature 
is outlined in red (dotted line). Bottom: Stars in fields Hll, fll6, 
and fl23 that fall within the triangular outline in the top panel 
are colored red. In the fields in which the ~ —300 km s~ x cold 
component is present, it is spread evenly as a function of projected 
distance from the minor axis. 

in §[731 

6.2. Metallicity Distribution 

So far we have considered only the kinematic properties 
of the minor axis population. The distribution of stellar 
metallicities, however, is also a powerful diagnostic of 
the presence and origin of substructure, since different 
galactic components (e.g., disk, inner spheroid) and tidal 
debris have different formation histories, and therefore 
different chemical abundances. 

Figure Qj}](a) displays [Fe/H] vs. w ne i for the M31 
RGB stars. The [Fe/H] values are based on a compar- 
ison of the star's position within the (I, V — I) CMD 
to a finely spaced grid of theoretical 12 Gyr, [q/Fe] = 
stella r isochrones (|Kalirai et al.ll2006bl : lVandenberg et alj 
2006) adj usted to the distance of M31 (783 kpc: Stanek & 
GarnavichTl^; iHohandl [199J). Stars with [Fe/H]< -1 
appear to be evenly distributed in velocity, while there is 
an obvious clump of metal-rich stars with velocities near 
—300 km s _1 . The bottom panels (b-d) show velocity 
histograms for stars in three [Fe/H] bins: (b) [Fe/H]> 
-0.5, (c) -1.0 <[Fe/H]< -0.5, and (d) [Fe/H]< -1.0. 
The strength of the cold component in each metallic- 
ity bin is measured by performing a maximum-likelihood 
double-Gaussian fit to the velocity distribution. Only 
the fraction of stars in the cold component (N su b/Nt t) 
is allowed to vary; the rest of the parameters are held 
fixed at the best-fit values from the fit to the full M31 
RGB sample (§I4.1| Table [2]). The fraction of stars in 
the cold component is 29.2% in (b), 20.9% in (c), and is 
negligible in (d), indicating that the cold component is 
metal-rich relative to the dynamically hot component. 

Figure [11] compares the metallicity distributions of 
stars in the velocity range of the ~ —300 km s _1 



Fig. 10. — (a) Metallicity vs. heliocentric velocity for the 
full M31 RGB sample. The majority of the population is metal- 
rich, with an evenly distributed metal-poor tail over the full range 
of velocities. A concentration of metal-rich stars near u^el = 
—300 km s _1 can be seen. The bottom three panels show ve- 
locity histograms for subsets of the data: (b) [Fe/H]> —0.5, (b) 

-1.0 <[Fe/H]< -0.5, and (d) [Fe/H]< -1.0. The 300 km s" 1 

cold component is more dominant in the metal-rich samples. A 
maximum-likelihood double-Gaussian fit was performed for each 
subset of the data, with all of the parameters except iV su b/iVtot 
held fixed at the best-fit values for the complete M31 RGB sample 
fij l4.ll Fig. [J] Table [2J. The cold component comprises a negli- 
gible fraction of the population in (d), 20.9% of the population 
in (c), and 29.2% of the population in (6), indicating that the 
~ —300 km s _1 cold component is relatively metal-rich. 

cold component (solid line) discovered in fields Hll, 
fl 16, and fl23 and stars that are identified with the 
hot spheroidal component (dashed line) in those fields. 
The hot spheroidal distribution (w oute r) is based on stars 
whose velocities are greater than ±2cr™ b away from 
(w) sub . This minimizes contamination of the spheroidal 
component by stars associated with the ~ —300 km s _1 
cold component. An additional constraint on the w oute r 
distribution is that only stars with w n ci < —150 km s _1 
are included in order to avoid residual MW dwarf star 
contaminants that lie in the range w no i > —150 km s _1 
(§ A. 2). Stars within the velocity range of the cold com- 
ponent can only statistically be identified as belonging 
to the hot spheroid or cold component, thus it is not 
possible to identify an uncontaminated sample of the 
~ —300 km s _1 cold component. Stars with velocities 
within ±lal uh of {v) sub are used for the Winner [Fe/H] dis- 
tribution, in order to maximize the number of members 
of the cold component while minimizing the contribution 
of spheroid stars. 

The [Fe/H] distribution of the Uouter sample has a 
peak at lower [Fe/H] values and a larger metal-poor 
tail than the [Fe/H] distribution of the Winner sample. 
The mean (median) metallicity of the w ou ter sample is 
<[Fe/H]) mcan = -0.72 (([Fe/H]) mcd = -0.63), while 
the mean (median) metallicity of the Winner sample is 
<[Fe/H]) m ean = -0.55 «[Fe/H]) med = -0.49). A KS 
test returns a probability of 0.7% that the two distribu- 
tions are drawn from the same parent distribution. The 
Winner distribution is highly contaminated by spheroid 
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Fig. 11.— The [Fe/H] distribution of the 

dinner (solid line) and 

foutor (dashed line) components in fields Hll, fll6, and fl23, in 
histogram (top) and cumulative (bottom) form. The ^i nne r (cold 
component) sample is defined to be stars with velocities within the 
range (v) suh - a s v uh < v hcl < (v) sub + cr= ub , where (v) suh and al ub 
represent the best-fit narrow Gaussian components in each field. 
This range was chosen to maximize the percentage of substructure 
stars compared to hot spheroid stars. The i) ou tor (hot spheroid) 
sample consists of stars with velocities Vy^ e \ < (v) Buh — 2c^ ub or 
(ti)™ b +2a Buh < tihcl < —150 km s" 1 . This minimizes contamina- 
tion from substructure stars and MW dwarf stars (§ A. 2) in the 
hot spheroidal [Fe/H] distribution. The «i nner sample is slightly 
more metal rich than the Pouter sample. Since the Winner sample 
is contaminated by spheroid stars, the true difference in [Fc/H] 
between the substructure and spheroid populations is somewhat 
greater than indicated in this plot. The cumulative [Fe/H] dis- 
tribution of stars in M31 ' s GSS (from fields a t 33 and 21 kpc; 
Guhathakurta et al. 2006; Kalirai ct al. 2006a) is plotted in the 
bottom panel for comparison (thin dotted line), and will be dis- 
cussed in § 17,11 

stars even within ±lcr^ ub of (w) sub (Fig. [6]); the estimated 
contamination of the Winner sample by spheroid stars is 
32.5%. A statistical subtraction of the u ou ter distribution 
(scaled by the contamination rate) from the Winner distri- 
bution yields a distribution with a mean (median) metal- 
licity of ([Fe/H]) me an = -0.52 «[Fe/H]) mcd = -0.45). 

An important effect on the measured [Fe/H] values is 
the assumed age (12 Gyr) of the stellar population. This 
assumption is wrong for at least field Hll, which has 
been shown by deep HST / ACS imaging to contain a wide 
spre ad of stellar ages, ranging from 6-13.5 Gvr ( Brown 
et al. 20031. l2006bl T Our data show a sharp cutoff in the 
[Fc/H] distribution at [Fe/H] — 0, with a metal-poor tail 

that extends out to [ Fe/H] 1.5 to -2.0 (Fig.fTTJ). The 

iBrown et all (|2006bh HST/ ACS data show a metal-rich 
cutoff in the [Fc/H] distribution at [Fc/H] — 0.3, with a 
metal-poo r tail which ex t ends to [Fe/H] — —1.5 to —2.0 
(Fig. 9 of IBrown et all ((2006b)). The star formation 
history derived from the HST/ ACS CMDs show that 
the intermediate age population (6-9 Gyr) is metal- rich 
([Fe/H]- 0), while the old (10-14 Gyr) stellar popu- 
lation is relatively metal-poor. Thus, our assumed age 
of 12 Gyr will underestimate metallicities for stars with 
[Fe/H] values near solar, but is appropriate for the more 
metal-poor stars ([Fc/H]< —0.5). Varying the age be- 
tween 6 and 14 Gyr introduces an » 0.3-0.4 dex spread 
in the [Fe/H] values derived from the isochrone fitting 
(|Kalirai et al.|[2006bl ); after accounting for this offset at 
the metal-rich end of the [Fe/H] distribution, our [Fe/H] 



distribution is consistent with the IBrown et al.l (|2006bl ) 
result. 

The intrin s ic spread in ages in the spheroid found by 
IBrown et al.l (2003, 2006b) can have two possible effects 
on our comparison of the [Fe/H] distributions of the Winner 
and fouter samples. If all the stars in these fields have 
a common spread in ages regardless of their kinematical 
properties, the error in the [Fe/H] measurement intro- 
duced by assuming a uniform age for the population will 
cause a shift in the actual [Fe/H] values, but the rela- 
tive difference between the Winner and w outC r populations 
will not be greatly affected. If the stars associated with 
the ~ —300 km s _1 cold component are systematically 
younger than the underlying smooth spheroid popula- 
tion, the measurement of [Fe/H] for the Winner sample 
derived from the 12 Gyr isochrones will be biased to- 
wards low metallicities, and thus the intrinsic difference 
in metallicity between the two populations will be greater 
than shown. 

We have also assumed that all stars are at the same 
line-of-sight distance as M31's center. This is a valid ap- 
proximation for the inner spheroid: at i? pr oj = 20 kpc, 
the spread in line-of-sight distances is expected to be 
about ±20 kpc (a spread in apparent magnitude of 
±0.05 dex), which corresponds to a spread in [Fe/H] 
of approximately ±0.03 dex. If the cold component is 
systematically more (less) distant than M31's spheroid 
(§ 17. 1[) . the intrinsic difference in metallicity between the 
populations will be slightly smaller (greater). 

Although fields fll5 and fl35 do not show clear evi- 
dence of substructure (§ 14. 2p , they are both within the 
minor axis range of the — —300 km s _1 cold component 
(Fig. 13 §EU). If the cold components in fields Hll, fll6, 
and fl23 have the same physical origin, it is reasonable to 
postulate that there may be substructure in fields fll5 
and fl35 that is not detected by the fits to the radial 
velocity distributions. The stars in fields fll5 and fl35 
that fall within the triangular region marked in the up- 
per panel of Figure [9] have mean (median) [Fe/H] values 
of -0.67 (-0.57) and -0.56 (-0.46), respectively. Stars 
from fields fll5 and fl35 that have velocities both out- 
side the triangular region and w ne i < —150 km s^ 1 have 
mean (median) [Fe/H] values of —0.61 (—0.53) and —0.64 
(—0.59), respectively. The [Fe/H] distributions of stars in 
fields fl35 and fl 15 that are within the triangular region 
of the ~ —300 km s _1 cold component are consistent 
with both the winner and w ou ter distributions. The differ- 
ence between the substructure and spheroid metallicity 
distributions is small, and is only statistically significant 
when the three fields contributing to the ~ —300 km s _1 
cold component are combined (into the Wi, m er and w outer 
samples). The number of stars within a restricted ve- 
locity range in any given field is too small to support a 
statistically significant comparison. 

6.3. Comparison to Previous Observations 

The Hll field has been presented in previous papers as 
a smooth spheroid field that is w ell-described by a si n- 
gle, kinematically h ot component (jBrown et al.ll2006aH bl: 
iKalirai et al.ll2006af ). There are two related factors that 
have caused this field to be reinterpreted as containing 
substructure. First, the data set presented in this paper 
represents an order of magnitude increase in the sample 
size of confirmed M31 RGB stars over previous spectro- 
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scopic samples published by our group in the SE minor 
axis region of the inner spheroid (i? P roj < 30). Only the 
Hll and aO fields have been previously published, and the 
data published in Hll contained only a fraction of the 
full data set for that field, based on preliminary reduc- 
tions (§ 12.2. 2j) . This paper presents the full Hll data set 
and combines it with new data from neighboring fields. 
Second, the triangular shape of the ~ —300 km s _1 cold 
component in a plot of velocity vs. position along the mi- 
nor axis implies that the debris has a relatively large ve- 
locity dispersion in the Hll field. This makes the debris 
harder to detect against the broad underlying spheroid 
than if it were kincmatically colder. The increase in the 
sample size of stars with recovered velocities in the Hll 
field, coupled with the context provided by the veloc- 
ity distributions in neighboring fields, has proved to be 
crucial in detecting the substructure in Hll. 

Previous observational studies have suggested the pos- 
sibilit y of substructure along M31's s outheastern minor 
axis. iReitzel fc Guhathakurtal (|2002fl found a dynami- 
cally cold grouping of 4 metal-rich M31 RGB stars (out 
of ~ 35) along a southeastern minor axis field at i? pro j ~ 
19 kpc. This hint of substructure was strengthened by 
subsequent observations at 7 and 11 kpc along the south- 
eastern minor axis, whi ch increased the tota l M31 RGB 
sample to ~ 1 00 stars ( Guhathakurtal 120021 ) . The star- 
count maps in iFerguson et all (|2002f ) (their Fig. 2) also 
show hints of a population of metal-rich stars along the 
southeastern minor axis, and deep HST/ACS imaging 
has discovered a significant interme diate-age population 
in field Hll (5151 iBrown et al.ll2003l ). With the spectro- 
scopic sample of > 1000 stars presented in this paper, we 
are able to confidently identify a cold component along 
the southeastern minor axis and characterize its proper- 
ties. 

7. PHYSICAL ORIGIN OF THE COLD 
COMPONENT 

The observed M31 RGB population along the south- 
eastern minor axis exhibits a spatially varying kinemat- 
ically cold component, which has a higher mean metal- 
licity than the underlying inner spheroid population. 
A cold component with these properties could be part 
of the outskirts of M31's disrupted disk or debris left 
by disrupted satellites. The models of F07 predict de- 
bris stripped from the progenitor of the GSS should be 
present in these fields. (It is also possible, of course, that 
the observed substructure is satellite debris unrelated to 
the GSS.) This section examines both the continuation 
of the GSS (§OJ and M31's disturbed disk (§OJ as 
possible physical origins of the cold component. 

7.1. Relation to the Giant Southern Stream 

7.1.1. Model of a Recent Interaction 

Debris in the form of coherent shells has been observed 
in many elliptical galaxies; these shells are believed to 
be formed by the tidal disru ption of a satel lite galaxy 
on a nea rly radial orbit (e.g.. lSchweizerlfl 980: Hcrnquist 
& Quinn 119881 ; iBarnes fc Hernquistl ll992). F07 presents 
the hypothesis that the Northeast and Western "shelves" 
observed in M31 (Fig.[TJ also see Fig. 1 of F07) are a sim- 
ilar phenomenon to these shell systems, and have been 
created by the disruption of the progenitor of the GSS. 
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Fig. 12. — Projected sky position (in M31 centric coordinates 
£ and n) of tidal debris in the F07 simulations of the merger of 
a dwarf galaxy with M31. Particles approaching their first peri- 
centric passage are part of the GSS (green). Particles approaching 
their second pericentric passage form the Northeast shelf (red), 
and particles approaching their third pericentric passage form the 
Western shelf (magenta) identified in F07. Particles in blue are ap- 
proaching their fourth pericentric passage, and form a faint shelf 
feature which is predicted to be most easily visible in the southeast. 
The position of our spectroscopic masks are also shown; fields fl23 
and fl35 straddle the edge of the Southeast shelf (Fig. [2j . The two 
ma sks at g = 0.3°, r\ = —1.5° are in field H13s, which is discussed 
in §[7X3] and §[8] 

Shells have coherent velocities, and display a distinctive 
triangu lar shape in the fi„ a vs. R v ^\ plane ( Merrifield & 
Kuiiken Tl998l . F07). In general, as flproj approaches the 
boundary of the shell the spread in velocities approaches 
zero, with the mean velocity at the tip of the triangle 
expected to be at the systemic velocity of the system. 

F07 present an TV-body simulation of an accreting 
dwarf satellite within M31's potential. The simulations 
use a static bulge+disk +halo model which i s based on 
the M31 mass models in lGeehan et al.l (J2006) combined 
with the observed stellar d ensity distribution in the halo 
(jGuhathakurta et"aTll2005l ). and assumes an isotropic ve- 
locity distribution in the outer halo. The simulated satel- 
lite's physical and orbital properties have been chosen to 
reproduce the observed propert ies of the GSS and N orth- 
east shelf, using the methods of iFardal et al.l ([2006) . The 
simulations show that the orbit which reproduces these 
features also reproduces a photometric feature identified 
in F07 as the 'Western shelf" and a n observed stream of 
count er-rotating planetary nebulae ([Merrett et al.l 120031 
12006T ). 

Figure [12] shows the projected sky positions (in M31- 
centric coordinates) of the satellite particles from the 
F07 simulation. The particles are color-coded by shell, 
or equivalently, by the number of pericentric passages 
they have made. Green particles represent particles ap- 
proaching their first pericentric passage; they correspond 
to the observed GSS. Red particles correspond to the ob- 
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served Northeast shelf; they are between their first and 
second pericentric passages. Magenta particles represent 
the Western shelf identified in F07, and are between their 
second and third pericentric passages. The blue particles 
are between their third and fourth pericentric passages, 
and represent the "Southeast shelf" predicted by F07. 
This last feature is predicted to extend out to a radius 
of 18 kpc and is expected to be very faint, as it con- 
sists of particles further forward in the continuation of 
the stream than the more visible Northeast and Western 
shelves. This feature actually covers ~ 180° in position 
angle on the east side of M31, although it is likely to 
only be visible in the southeast due to its overlap with 
the Northeast shelf and M31's disk. In the F07 simula- 
tions, the Northeast shelf is made up of both the leading 
material from the progenitor's first pericentric passage 
and trailing material from its second pericentric passage, 
while the Western shelf is formed by leading material. 
The simulations are unable to constrain whether or not 
the satellite disrupts completely, as this is dependent on 
the central density of the satellite. 

The bottom panel of Figure [13] presents the distribu- 
tion of particles from the F07 simulation in the t>i os vs. 
-Rproj plane. The figure shows particles related to the 
merging satellite as well as particles associated with the 
static bulge+disk+stellar-halo M31 model used in F07. 
In order to carry out a precise comparison to our ob- 
servational data set, the F07 simulation particles were 
selected based on their projected sky position; all parti- 
cles that fall inside a 16'xlO' area (the approximate area 
of one DEIMOS mask is 16'x4') centered on the posi- 
tion of our observed fields, and oriented at the position 
angle of our observed masks, are displayed. The satel- 
lite particles are color-coded by shell (or, equivalcntly, 
the number of orbits they have made) as in Figure 1121 
Green particles are associated with the GSS, red parti- 
cles with the Northeast shelf, and black particles with 
the bulge+disk+stellar-halo model for M31. The blue 
particles, which are part of the predicted "Southeast" 
shelf, form the distinctive triangular shape expected of a 
shell feature in the -R pro j-vios plane. 

In the following discussion, it is important to distin- 
guish between "spillover" from the GSS (material asso- 
ciated with the GSS that is in roughly the same orbital 
phase as the material in the GSS, i.e., green particles in 
Figs. [12] & IT3|) versus wrapped around portions of the 
GSS (material associated with the GSS progenitor that 
is leading the GSS and has undergone one or more ad- 
ditional pericentric passages, i.e., red, magenta and blue 
particles in Figs. I12lfe ll3[) . The latter is the main theme 
of this paper, although we also briefly discuss the former 
in §§7.1.2-7.1.3. 

7.1.2. Comparison of Model to Data 

SKY POSITION 

Figure [12] shows the projected sky positions (in M31- 
ccntric coordinates) of the satellite particles from the F07 
simulation as well as the size, position and orientation of 
our Keck/DEIMOS slitmasks (rectangles). Fields fl23 
and fl35 land on the edge of the predicted Southeast 
shelf (blue particles), fields fll5, fll6, Hll and fl09 all 
lie within the boundary of the Southeast shelf, and fields 
fl30 and aO lie beyond it. Thus, the model predicts that 
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Fig. 13. — Line-of-sight velocity vs. projected radial distance 
from M31's center (R pI0 j) of spectroscopically confirmed M31 RGB 
stars (top panel) and particles from the F07 simulations of the orbit 
of the progenitor of the GSS (bottom panel). Particles are drawn 
from the locations of each of our DEIMOS masks. The area from 
which particles have been drawn has been increased relative to the 
size of a slitmask to increase the number of particles. The satellite 
particles are color-coded according to which shell they are in: the 
giant southern stream (green), the Northeast shelf (red), and the 
predicted Southeast shelf (blue). Black points are particles from 
the bulge+disk+stellar-halo of M31. The blue particles form a 
triangular shape, with an increasingly wide kinematic profile as 
the minor axis distance to the center of M31 decreases, as seen 



in the data. The tip of the triangle at R p 



1.3° (18 kpc) in 



the simulated data agrees well with the observed tip in the data 
in field fl23. The cold concentration of M31 particles at i>M3l 



km s 1 extending from R pv , 



< 0.58° to ~ 0.73° (7.9 to 10.0 kpc) 



corresponds to the disk of M31 (bottom panel); the kinematical 
signature of a smooth, cold disk is not seen in our data (top panel). 

the edge of the shell feature should pass directly through 
our CFHT/MegaCam image. Indeed, there is an ap- 
parent edge visible in the CFHT starcount map (Fig. [2] 
passing through field fl23), in the same location as that 
pre dicted for the South east shelf. A close inspection of 
the llbata et alj (|2005f ) starcount map (Fig. [1]) reveals 
a point of bifurcation between the edge of the North- 
east shelf and a fainter feature at £ « 1.6°, rj « 0.2°, in 
rough agreement with the bifurcation of the two features 
in Figure 1121 This bifurcation is more evident in the 
Sobel-filtered map in F07 (their Fig. 1). The radii of the 
shells in the simulation are robust (§ 4.2 of F07); thus the 
agreement between the observations and the simulations 
is a strong confirmation of the validity of the F07 model. 

KINEMATIC TRENDS 

The ~ —300 km s" 1 cold component observed in our 
minor axis fields shows the distinctive triangular veloc- 
ity pattern expected of a shell feature in the i? pr oj-vios 
plane. Figure IT31 compares our data (top panel) to the 
F07 model (bottom panel). The distribution of observed 
velocities narrows to a tip at i? pro j ~ —1-3° (18 kpc) in 
the simulated particle distribution, which is similar to 
the position of the tip of the velocity distribution in our 
observed data. At i? pro j ~ — 1° (13.7 kpc), the velocity 
distribution of the observed substructure has widened to 
a spread of ~ 200 km s" 1 (measured from the edges of 
the feature), also in agreement with the velocity spread 
of the predicted Southeast shelf. 

The "boxy" shape of the velocity distribution in field 
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Hll (Fig. [6]) is also consistent with the interpretation 
of the substructure as being part of a shell system. The 
velocity distributions of shells have a clearly defined min- 
imum and maximum line-of-sight velocity at a given 
i?p r oj. Stars tend to congregate at the minimu m and 
maximum velocities (|Merrifield fc Kuijkenl fl998l . F07), 
although their location in the i? pr oj-^ios plane depends 
on the region they occupy in space (cf. the discussion in 
F07). 

A maximum-likelihood Gaussian fit to the particles 
identified with the Southeast shelf and within the mi- 
nor axis distance spanned by field fl23 yields parameters 
of (v) = -280.6 km s^ 1 and a v = 19.4 km s -1 . A 
Gaussian fit to the Southeast shelf particles within the 
minor axis range spanned by fields fl 16 and Hll returns 
(v) = —292.6 km s _1 and a v — 60.5 km s _1 . The mean 
velocity and dispersion of the predicted shelf is in good 
agreement with the properties of the observed substruc- 
ture (Table [2]). 

METALLICITY DISTRIBUTION 

If the substructure identified in this paper is part of the 
predicted Southeast shelf in F07, it should have a similar 
metallicity distribution to that of the GSS, since the two 
structures originated from the same progenitor. As part 
of our Keck/DEIMOS survey of M31's inner spheroid 
and halo, we have taken spectra in two fields located on 
the G SS: a field at -R pro j = 33 kpc (|Guha thakurt a~et al.1 
2006) and a field at i? proj = 21 kpc (H13s: iKalirai et all 
2006a). The cumulative [Fe/H] distribution of stars iden- 
tified kinematically as belonging to the GSS in these two 
fields is plotted in the bottom panel of Figure [IT] (thin 
dotted line). It is very similar to the [Fe/H] distribution 
of stars that are kinematically associated with the sub- 
structure in fields Hll, fll6, and fl23. The mean and 
median [Fe/H] of the stars in the GSS are 0.1 and 0.05 
dex more metal-poor, respectively, than the mean and 
median [Fe/H] of the substructure in fields Hll, fll6, and 
fl23, after correcting for spheroid contamination (§ 16. 2|) . 
The estimated number of inner spheroid star contami- 
nants iin_JbeGtSSj3amr^^ 

et al. 120061 : IKalirai et alj T2006a). If the 3 most metal- 
poor stars are removed ([Fe/H]< —2.25) from the GSS 
distribution, the mean and median metallicity of the GSS 
stars are only 0.01 dex more metal-poor than the the mi- 
nor axis substructure. The [Fe/H] values of the GSS 
stars have not been corrected for the GSS' measured dis- 
tance relative to M31 (~ 50 kpc behind M31 for these 2 
fields: iMcConnachie et al.ll2003l ). Accounting for this ef- 
fect would decrease the average metallicity of the GSS by 
~ 0.1 dex. The distance to the minor axis substructure 
is not known, although the F07 simulations predict that 
the Southeast shelf should be approximately at M31's 
distance, with a spread in distances of ±9.2 kpc (this 
corresponds to ±2<r in terms of the distribution of parti- 
cle distances). 

STRENGTH OF THE COLD COMPONENT 

The cold component comprises 44% of the total popu- 
lation of observed stars in fields Hll and fll6 and 31% 
of observed M31 RGB stars in field fl23 (Table©. This 
corresponds to a lower limit for the total fraction of stars 
in the cold component of 21.7% in the fields within the 
predicted range of the Southeast shelf (fl09, Hll, fll6, 



fll5, fl23, and fl35). The Southeast shelf in the sim- 
ulations is much weaker, comprising only 2.7% of the 
total population in these fields (this number increases to 
3.4% if the number of shelf particles is compared only to 
the number of bulge+disk+stellar-halo M31 particles). 
The strength of the feature in the simulations is highly 
dependent on the mass of the progenitor and the time 
since the first collision (F07). Thus, the strength of the 
observed substructure will place interesting constraints 
on future models of the stream, but cannot be used as a 
reliable discriminant of the applicability of the model at 
the present time. 

FIELDS WITHOUT CLEAR DETECTION OF 
SUBSTRUCTURE 

We do not find a clear detection of substructure in 
fields fl30 and aO. In the context of the Southeast shelf, 
this is not surprising as both these fields are beyond the 
radial range spanned by the shelf (Fig. [T2|) . Field fl09 
is significantly inward of the innermost field in which 
we detect substructure. In the simulation, the particles 
associated with the Southeast shelf continue into the re- 
gion covered by field fl09 with a spread in velocities of 
~ 350 km s" 1 (Fig. [131 bottom panel). The data in 
field fl09 is consistent with being drawn from a single 
Gaussian (§ I4.2|I and shows no evidence of substructure. 
A secondary component with a spread in velocities as 
wide as predicted would be very difficult to differentiate 
from the broad spheroidal component, and would require 
a much larger sample of M31 RGB stars in this field than 
is currently available. 

Field fl 15 is well within the boundaries of the South- 
east shelf. As discussed in § 16. 1( its velocity distribution 
is consistent with being drawn from the same parent dis- 
tribution as field fl 16. The shell in field fl 15 may be diffi- 
cult to detect in our data due to the broad (~ 55 km s _1 ) 
nature of the substructure and the smaller number of 
stars available in this field (~ 50% less than in field fll6), 
or the shelf may be inherently clumpy. 

Field fl35 is on the edge of the simulated shelf, 
and shows evidence for a peak of stars near Uhd ~ 
—300 km s _1 in its radial velocity histogram. In light of 
the simulations, we discuss this field in detail in § I7.1.3| 
and show that it may contain a kinematically-cold com- 
ponent whose properties are consistent with both the ob- 
servations of the ~ —300 km s~ x cold component in fields 
Hll, fll6, and fl23 and the simulations of the Southeast 
shelf. 

DEBRIS FROM THE GSS 

The GSS is observed to have an asymmetric shape, 
with a sharp edge on the eastern side and a more grad- 
ual d ecline in d ensity o n the western side f McConnachie 
et al. I2003D . However, the eastern edge is not an ab- 
solute one in the models, and "spillover" material from 
the GSS (§ 17. 1[) is predicted by the F07 simulations to be 
present in all of our fields. Field fl35 is the closest of our 
fields to the GSS' eastern edge (Figs, ffl and [12]); § 17X31 
discusses the evidence for spillover debris from the GSS 
in this field. 

Since the density of M31 RGB spheroid stars falls off 
strongly with increasing radius in the inner spheroid, the 
contrast of cold GSS debris against the dynamically hot 
spheroid is expected to be greatest in our outermost field, 
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aO. The GSS debris in field aO is predicted to have a 
mean velocity of (v)gss — —364 km s _1 . Although field 
aO shows hints of peaks in the radial velocity distribution 
at Uhci < —300 km s _1 (Fig. [6]), only a handful of these 
stars are as metal-rich as the GSS ([Fe/H]> —0.75). This 
allows us to place an upper limit on the contamination 
of field aO by the GSS of < 5 stars (< 6%). 

Although the F07 simulations reproduce many of the 
observations in the GSS, they predict a much larger 
amount of debris on the eastern side than is observed 
in our fields (Fig. [T3|) . Many factors can influence the 
structure of the debris in the simulations, including the 
shape and rotation of the progenitor. The current mod- 
els of the stream (F07) use a spherical, non-rotating pro- 
genitor. A more complex model of the progenitor may 
be required to reproduce the observations (Fardal et al., 
private communication) . 

7.1.3. Evidence of the Southeast Shelf and Spillover Debris 
from the GSS in Field fl35 

The F07 simulations predict that in addition to the 
Southeast shelf some spillover debris associated with the 
GSS (§[71]) should be present in field fl35. The sky co- 
ordinates of field fl35 are £ = 0.7, r\ = —1.1, which 
places it on the edge of both the GSS and the Southeast 
shelf in the F07 simulations ( Figure [T2"j) . In addition, 
the radial velocity distribution of stars in field fl35 is 
not well-fit by either a single or double Gaussian f? l4.2[) 
and shows evidence of a metal-rich population (§ 16. 2p . 
Motivated by the close match between the observations 
and simulations of substructure in fields Hll, fll6 and 
fl23, we carry out a constrained fit of the radial velocity 
histogram of field fl35 to the sum of three Gaussians to 
determine if the Southeast shelf is present in this field. 

The mean velocity and velocity dispersion of the sim- 
ulated Southeast shelf and GSS particles in field fl35 
are (v) SE = -286 km s~\ cr^ E = 19 km s" 1 and 
(w)gss = -458 km s -1 , cr° ss = 40 km s" 1 . These val- 
ues were used as rough constraints for the triple Gaussian 
fit: the means were allowed to vary within ±100 km s _1 
of the predicted values and the dispersions were al- 
lowed to vary from 1 km s" 1 to 3er^ E and 2<r^ ss . The 
wide Gaussian component parameters were held fixed 
at the values for G sph (v): (v) sph = -287.2 km s' 1 and 
<P h = 128.9 km s" 1 (§EQ). The maximum-likelihood 
triple Gaussian fit is displayed in Figure fMl (solid curve). 
The wide underlying inner spheroid component (G sph (v), 
dot-dashed curve) comprises 45% of the population. The 
Southeast shelf component (dashed curve), which is the 
narrow peak at (v) sub = —273 km s _1 , comprises 30% of 
the population and has a width of cx^ ub = 30 km s _1 . The 
"GSS" component (dotted curve) at (v)= —449 km s _1 
has a dispersion of a v — 55 km s~ x and comprises 25% 
of the total population. If a more constrained fit is car- 
ried out with the (v) and o~ v parameters for all three 
Gaussian components held fixed (at the predicted values 
for the simulated shelf and stream particles and at the 
parameters of the Gaussain G sph (?;)) and only the frac- 
tions of stars in the various components are allowed to 
vary, the best-fit distribution has iV s heif/-/Vtot = 0.18 and 
Noss/Ntot = 0.17. 

The kinematic properties (mean velocity and velocity 
dispersion) of the Southeast shelf component in the triple 
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Fig. 14. — Radial velocity histogram of M31 RGB stars in field 
fl35. A constrained triple Gaussian (solid curve) has been fit to the 
observed data using a maximum-likelihood technique, with rough 
constraints imposed on the p arame ters based on the properties of 
the simulated substructure (§ 17.1.30 . The observed velocity distri- 
bution is well fit by a sum of three Gaussians: (i) G Bph (v), the 
wide Gaussian which corresponds to the underlying inner spheroid 
of M31 (dot-dashed curve, § 14.11 1. (ii) a component centered at 
(v) suh = -273 km s" 1 with a width of o-J ub = 30 km s"\ which 
comprises 30% of the total population and which likely corresponds 
to the Southeast shelf (dashed curve), and (iii) a narrow component 
centered at (v)= —449 km s — 1 with a width of a v = 55 km s _1 , 
which comprises 25% of the total population (dotted curve). The 
mean velocity and velocity dispersion (±l<x„) of the cold compo- 
nents in field H13s, at a similar radial distance along the GSS as 
field fl35, arc shown as arrows and horizontal lines (§17.1.31 Kalirai 
et al. l2006al) . 



Gaussian fit to field fl35 are consistent not only with the 
simulations, but also with what one would expect for 
the Southeast shelf in this field based on the observa- 
tions (e.g., Fig. [13l Table [2]). In addition, the fraction 
of the population in fl35 which is in this component is 
consistent with the fraction of the population which is 
identified with the Southeast shelf in fields Hll, fll6, 
and fl23. As further evidence that the Southeast shelf 
is detected in field fl35, Figure 15 shows [Fe/H] vs. <>i ic i 
for the M31 RGB stars in field fl35 (panel a), as well 
as velocity histograms for stars in two [Fe/H] bins [(b) 
[Fe/H]> -0.75 and (c) [Fe/H]< -0.75]. As in fields Hll, 
fll6 and fl23, the substructure in field fl35 that is iden- 
tified with the Southeast shelf ((v}= —273km s _1 ) is 
relatively metal-rich. 

The velocity dispersion of a v = 55 km s _1 for the "GSS 
component" inferred from the first of the triple-Gaussian 
fits above is large compared to previous meas urements of 
~ 15 km s" 1 for the dispersion of the GSS (llbata et al.l 
120041 : iGuhathakurta et al.ll2006l : iKalirai et al.ll2006al ). In 
addition, although the best-fit {v) suh of the most negative 
component in fl35 is similar to that of the F07 model, the 
predicted GSS mean velocities in the F07 model are not 
as ne gative as the observed velocities of the GSS. Kalirai 
et al. (|2006af ) analyzed a field centered on a high sur- 
face brightness protion of the GSS, at approximately the 
same radial distance along the stream as field fl35 (field 
H13s, located at £ = 0.3°, 77 = -1.5° in Figure [TJJ) and 
found a secondary cold component, the "H13s secondary 
stream," whose origin and physical extent are unknown. 
The observed mean velocities of the GSS and secondary 
stream in the H13s field are (v) GSS — —513 km s -1 and 
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Fig. 15. — (a) Metallicity vs. heliocentric velocity for the 
M31 RGB stars in field fl35, which lies to the east of the edge of 
the GSS. The mean velocities and dispersions (±la v ) of the GSS 
and th e secondary stream from the GSS field H13s IIKalirai et al.l 
2006a) are marked, as in Figure [14] (&) Velocity distribution of 
stars with [Fe/H]> -0.75. In addition to the ~ -300 km s" 1 cold 
component, the metal-rich subset shows evidence of concentrations 
of stars corresponding to the GSS and secondary stream in the H13s 
field, (c) Velocity distribution of stars with [Fe/H]< -0.75. The 
velocities of stars with metallicities lower than —0.75 dex appear 
evenly distributed. 

(u) scc str - = —417 km s -1 , respectively, with velocity dis- 
persions of cr^ ss = cr^ cc str - = 16 km s" 1 (Kalir aTet al.1 
20061). The mean velocities and velocity dispersions 
(±lcr„) of these two components are shown as arrows 
with horizontal lines in Figures [HI and [T5l There appear 
to be two metal-rich peaks in fl35 with the approximate 
velocities of the GSS and secondary stream in H13s. If 
both the GSS and H13s secondary stream are present 
in field fl35, they appear in approximately equal pro- 
portion. In H13s the GSS dominates over the s econdary 
stream by a factor of two dKalirai et al.|[2006af ). 

In the simulations, the GSS dominates over the South- 
east shelf in field fl35 by about a factor of 10; the ob- 
servations indicate that, at best, these two populations 
are roughly equal. The simulations predict that the GSS 
should comprise a total of 35% of the population in field 
fl35, which is somewhat larger than the fraction of the 
total population in the most negative cold component 
from the triple-Gaussian fit shown in Figure [TJ] (25%). If 
the "GSS" component from the triple-Gaussian fit is ac- 
tually comprised of two narrower streams, the GSS com- 
prises a much smaller percentage of the stars in this field 
than predicted by the simulations. 

7.2. Arguments Against a Disk Origin 

The substructure discovered in fields Hll, fl 16 and 
fl23 is centered at close to the systemic velocity of M31 
(u sys = —300 km s^ 1 ), which is also the radial velocity 
expected for an M31 disk component on the minor axis. 
Recent observational evidence suggests that M31's stellar 
disk extends smoothly out to i?disk ~ 40 kpc and has 
a velocity dispersion of ~ 30 km s _1 ; isolated features 
with disk-like kine matics have been observed as far ou t 
as -R disk ~ 70 kpc (|Reitzel et al.ll2004l ; Hbata et al.ll2005l ). 



We present three lines of evidence arguing against a disk 
origin for the ~ —300 km s" 1 cold component found on 
the southeastern minor axis. 

DISK TO INNER SPHEROID SURFACE BRIGHTNESS 
RATIO 

Previou s measurements of M31's stella r disk and inner 
spheroid (IWalterbos fc Kennicuttlll988l; Pritchet & van 
den Bergh|1994j) indicate that the disk constitutes a neg- 
ligible fraction of the total light in our minor axis fields. 
Even in our innermost field fl09 (at i? pro j ~ 9 kpc, cor- 
responding to i?disk ~ 38 kpc for a disk inclination of 
77°), the disk fraction is expected to be only ~ 10% 
(|Guhathakurta et al.ll2005l) base d on disk scale radii of 
5.0-6.0 kpc (jlbata et aTTl2005f ). We see no evidence 
of a cold, disk-like feature in fl09's radial velocity dis- 
tribution, although a < 10% component could be diffi- 
cult to detect. The disk fraction drops sharply at larger 
radii: for example, at the distance of our innermost field 
containing the ~ —300 km s _1 cold component, Hll 
(-Rdisk ~ 51 kpc), the exp ected smooth disk fraction is 
~ 1% dBrown et al.ll2006ah . 

Non-uniformities in M31's stellar disk could result in a 
higher disk fraction in our fields. To explain the strength 



of the 



-300 km s cold component in field Hll 



(N su t)/N to t — 44%), the disk would have to contain a 
45 x enhancement in this field relative to the smooth 
disk. Even more extreme disk enhancements are needed 
to explain the cold component in fields fll6 and fl23. 

If there is a warp suc h that the outer di sk is more face 
on than the inner disk (Iba ta et al.l {2005) find a best fit 
disk inclination angle of 64.7° from 20 < -Rdisk < 40 kpc), 
the effective disk radii of our fields will be smaller and 
the smooth disk contribution larger than in the above 
calculation. As an extreme example of a warp, we con- 
sider a disk of scale length 5.7 kpc whose inclination 
changes from 77° at small radii ( -Rdisk < 20 kpc) to 60° 
at large radii (-Rdisk > 20 kpc). In this case, Hll is at 
-Rdisk ~ 35 kpc while fl23 is at i?di S k ~ 47 kpc. With 
this disk model, taking into account projection effects 
in the disk surface brightness (the surface brightness de- 
creases as the disk becomes more face on) but ignoring 
dust effects, the expected disk fractions in fields fl09, 
Hll, fll6, and £123 are 20%, 18%, 17%, and 10%, re- 
spectively. Thus, it would still require a ~ 2.5 x enhance- 
ment in the disk to explain the cold component in Hll. 
More importantly, the 20% cold, smooth disk fraction 
predicted by this warp model in field fl09 is inconsistent 
with our radial velocity data (Fig. [5]). This disk fraction 
estimate for the warp is optimistic (high) in that it does 
not account for the reduction in surface brightness that 
would be caused by any stretching associated with the 
putative warp. 

VELOCITY DISPERSION 

The measured velocity dispersion of the 
~ —300 km s" 1 cold component in fields Hll and 
fl 16 is > 50 km s _1 , while the velocity dispersion 
of this component is only ~ 10 km s _1 in field fl23. 
These measurements are significantly above and below, 
respectively, the typical velocity dispersion of 30 km s _1 
meas ured for the ex tended, disk-like structure by Ibata 
et al. (|2005f) . for which thev observe a range in velocity 
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dispersions from ~ 20-40 km s _1 . If the cold com- 
ponent in each of these three fields is from the disk, it 
would require a warp and non-uniformities in the disk 
that happen to have the observed triangular shape in 
the position- velocity plane (top panel of Figure [9|) . 

STELLAR AGES AND METALLICITIES 

If the ~ —300 km s _1 cold component is debris as- 
sociated with M31's extended disk, its stellar popula- 
tion should reflect this. However, a comparison of the 
stellar ages and metallicities of the Hll field with deep 
HST/ACS and ground-based imaging of disk-dominated 
M31 fields yields very different star formation histories. 
Obviously disturbed sections of M31's disk (e.g., the 
Northern Spur and the Gl clump) show evidence of re- 
cent star formation (~ 3 Gyr ago in th e Northern Spur 
and ~ 250 Myr ago in the Gl clump; iFerguson et al.l 
US), while deep HST/ACS imaging of field Hll re- 
veals that very few of the s t ars are younge r than 4 Gyr 
(IBrown et all [20031 . I2006bl ). IBrown et all (l2006bl ) find 
that their disk-dominated HST/ACS field H13d, located 
at 25 kpc along the northeastern major axis, contains 
a significantly younger and more metal-rich p opulation 
than Hll. In addition, IBrown et alJ (|2006al ) find very 
good agreement between the stellar populations of field 
Hll and an HST/ACS field on the CSS (H13s, discussed 
further in §[8]). 

The comparison of the stellar ages and metallicities of 
M31 disk fields vs. that of field Hll is complicated by 
several uncertainties. A radial gradient in disk prop- 
erties could result in M31's outer disk (in field Hll) 
being more metal-poo r and older than the inner parts 
of the disk. Also, the IBrown et~"aTI (|2006bD HST/ACS 
field H13d likely includes multiple galactic components 
including spheroid and wrap-around debris from the pro- 
genitor of the GSS, although the radial velocity distribu- 
tion indicates that M31's disk is the dominant component 
in this field (see Fig. 9 of IKalirai et~aTll2006al ). 

In conclusion, the observations disfavor an extended 
rotating disk model as the physical origin of the ~ 
—300 km s _1 cold component identified along M31's 
southeastern minor axis. Although a disk origin for this 
substructure cannot be ruled out, it requires simulta- 
neous contrivance of multiple properties of the disk to 
explain the observations: patchy disk structure (to ex- 
plain the absence of an observed disk in fl09), a warp 
and enhancement (which cannot be due to recent star 
formation) of the disk to explain the strength of the cold 
component, an anomalously large velocity dispersion for 
fields Hll and fll6 and a smaller than average velocity 
dispersion in field fl23, and a significant radial gradient 
in the metallicity and age of the stellar disk populations. 
Compared to the elegance of the southeastern shelf in- 
terpretation based on F07's simulations, which is a true 
prediction and explains the observed properties of the 
cold component in all the fields in which it is detected, 
the disk origin clearly fails the test of Occam's razor for 
the most likely physical origin of the ~ —300 km s _1 
cold component. 

8. IMPLICATIONS FOR THE 
INTERMEDIATE- AGE SPHEROID 
POPULATION 
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(f2003l . [2006al lbh present HST/ACS pho- 
tometry of fields in M31 down to 1 — 1.5 magnitudes be- 
low the main- sequ ence turnoff. Our f ield Hll is co- 
incident with the IBrown et al.l (|2003f ) spheroid field 
Fig. [2 l). Tha photometry presented in IBrown et al.l 
2006al ) is from a field on the GSS of Andromeda at a 
projected radial distance of 20 kpc, and is coincident with 
the Keck/DEIMOS s pectroscopy field H13s presented in 
IKalirai et alJ (|ji006al ) and shown in Figure [T2l = 0.3°. 
■ n = — 1.5°). In the "smooth" spheroid field. IBrown et aTl 
((2003h find th it - 30% (by mass) of the stellar pop- 
ulation is intermediate-age (6-8 Gyr) and metal-rich, 
while another 30% o f the population is old (11-13.5 
Gyr) and metal-poor. IBrown et al. (2006a]) find remark- 
able agreement in the CMDs of the stream and spheroid 
fields, indicating that the two fields have very similar age 
and metallicity distributions. They query whether the 
similarities bet ween the populations could "be explained 
passing through the spheroid field," but 
explanation is problematic: the stream 
would have to dominate the spheroid by approximately 
the same factoi in both fields (3:1 based on a kinematica l 
stream field H13s: IKalirai et "aTl 12006a). 
ical profiles of the two fields are distinctly 
the Hll field failing to show the cold 
i -1 ) signature of the stream seen in the 
H13s field. However, they presciently suggested that the 
similarity in tl e two populations (spheroid and stream) 
implies that " ;he inner spheroid is largely polluted by 
material stripped from cither the stream's progenitor or 
similar objects " 

In light of t re substructure presented in this paper, 
this seems to bs the correct interpretation of the similar- 
ity between the "spheroid" and GSS stellar populations. 
The spatial and kinematic properties of the substructure 
suggest that the r egion of the spher oid imaged in the 
orig inal HST ffeld {B rown et al. 2003) is in fact contami- 
nated by stars from the progenitor of the GSS. The kine- 
matical signature of the substructure at the minor axis 
distance of the HST/ACS field (Hll) is both predicted 
and observed ',o be relatively wide (Fig. [T3j) . and thus 
less obvious against the underlying hot component. In 
the context of the F07 simulations, the minor axis sub- 
isolated, but is part of one of a series of 
y the disruption of the GSS' progenitor, 
which collectively contaminate a large part of the inner 
spheroid of MS 1 (Fig. [TJJ). 
The current analysis suggests that ~ 45% of the M31 
the Hll field are in fact part of the 
1 cold component, and not part of the 
broad spheroid. In the H13s field, 75% of the M31 RGB 
stars are part of a cold component (|Kalirai et a l. 2006a) . 
This difference in substructure fraction agrees nicely with 
the difference n the fraction of intermediate-age (< 10 
Gyr), metal-rich stars found in the stream and spheroid 
fields in lRrowIi et all ([2006b): 70% vs. 40%, respectively. 
However, recent HST/ACS observations of a field in the 
location of oui fl30 mask s at 21 kpc imply t hat this is 
not the end of the story: IBrown et al.l (|2007l ) find that 
the stellar population in Hll can not be fit by a linear 
combination oi the GSS (H13s) and the 21 kpc spheroid 
(fl30) stellar populations, due largely to the presence of a 
greater number of stars younger than 8 Gyr in Hll than 
in the GSS fiqld. Nevertheless, the observational cvi- 
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dence, combined with the theoretical predictions of F07, 
strongly favor the explanation that the age and metallic- 
ity distributions of the stream and spheroid HST fields 
are so remarkably similar because the same progenitor 
polluted both fields with substructure. 

9. SUMMARY 

The use of the diagnostic method described in Gilbert 
et al. (|2006f ) has enabled us to isolate the first sample 
of spectroscopically confirmed M31 RGB stars defined 
without the use of radial velocity. We use this sample 
of ~ 1000 M31 RGB stars to measure the velocity dis- 
persion of the inner spheroid of M31; in the radial range 
-Rproj =9 — 30 kpc the inner spheroid has a velocity dis- 
persion of er^ ph =129 km s _1 . Our data show no evidence 
of a decrease in the velocity dispersion over this radial 
range. 

The stellar radial velocity distribution in these fields 
shows evidence of a significant amount of substructure. 
Compared to the large velocity dispersion seen in the 
underlying hot spheroid population, the substructure is 
kincmatically cold, exhibiting a decrease in velocity dis- 
persion with increasing projected radius. In the fields in 
which the ~ —300 km s _1 cold component is observed, 
w 41% of the stars are estimated to belong to it; the 
rest are members of the hot inner spheroid of M31. The 
metallicity of the substructure is higher than that of the 
broad spheroidal component in the fields in which it is 
observed. 

The physical origin of the substructure discovered in 
this paper is most likely tidal debris stripped from the 
progenitor of the GSS. The data agree very well with 
the location and kinematical properties of the Southeast 
shelf predicted by the F07 simulations of the disruption 
of the GSS' progenitor, and will add significant obser- 
vational constraints to those already existing from the 
GSS, Northeast shelf, and Western shelf, enabling de- 
tailed modeling of M31's dark matter distribution (F07). 
The minor axis fields also place constraints on the spatial 
distribution of the GSS itself. The GSS contamination in 
our minor axis fields is much smaller than predicted by 
the current models of the stream (F07), which suggests 
the stream's progenitor had a more complex structure 
than the spherical, non-rotating models used so far. 

The newly-discovered substructure sheds light on the 
discovery of a significant interme diate-age popul a tion i n 
the "smooth" spheroid field by iBrown et al.l (2003), 



and the subsequent discovery of the similarity in ages 
and metallicities of the stars in the sp heroid field and a 
field on the GSS ([Brown et alJ l2006al fbh . The spheroid 
HST/ ACS field was not in fact placed on a "smooth" 
spheroid field, and the intermediate-age population may 
be part of the substructure observed in this field. If the 
substructure identified in this paper is indeed from the 
same progenitor as the giant southern stream, it is not 
surprising that the two HST/ ACS fields would have very 
similar age and metallicity distributions. Given the num- 
ber of observed fields in the inner spheroid which are 
contaminated by subs tructure, both in the current work 
and in the literature (I .Irwin et all 120051 : iFerguson et~all 
120051 ; iKalirai et alj|2006al ), it seems likely that the in- 
ner spheroid is highly contaminated by tidal debris. A 
"smooth" inner spheroid field may in fact be a rarity. 
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APPENDIX 

These appendices are intended to give the interested reader more insight into the origin of the differences in the 
distributions of the M31 RGB samples shown in Figure [3] (§ A.l). The amount of contamination vs. completeness for 
different (Li) thresholds is discussed in § A.2. Finally, in § A. 3 we quantify the effect of the dwarf contamination on 
the radial velocity distributions for the M31 RGB sample used in this analysis. 

A.l Effect of the Radial Velocity Diagnostic 

As discussed in §[3j empirical probability distribution functions (PDFs) based on training sets of M31 RGB and MW 
dwarf stars are used to determine the probability that an individual star is an M31 red giant (-Pgiant) or MW dwarf 
(Pdwarf) in 4 (5 for field aO) diagnostics. These probabilities are used to compute the likelihood a star i is a red giant 
in a given diagnostic j: 



(1) 
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Fig. 16. — Overall likelihood distributions with ((Li), dashed) and without ((Lj)y, solid) the inclusion of the radial velocity diagnostic. 
The analysis in this paper is based on the sample of M31 RGB stars shown by the shaded histograms. For stars with {£j)> 0, removal 
of the radial velocity diagnostic from the (Li) calculation causes a shift in the peak of the distribution from {Li)~ 2 to {Li)~ 1. This is 
expected, since radial velocity is the most powerful of the five diagnostics. A similar shift of the peak of the distribution towards (Li)= 
likely occurs for stars with (Li)< 0, however the number of stars in this category is too small for this shift to be easily discernible. 



A star's overall likelihood of being an M31 RGB star is defined as 

Yl w 3 L ij 

(Li) = -4= . (2) 

3 

All available diagnostics for a star are given a weight of one, unless the star is an outlier with respect to both the M31 
RGB a nd MW dwarf ' PDF s in a two-dimensional diagnostic, in which case the weight of that diagnostic is reduced 
(§3.3 of lGilbert et"aIIP2006h . 

Figure [TBI shows the overall likelihood distributions for each field, both with ((Li), dashed histogram) and without 
((Lj)y, solid histogram) the radial velocity diagnostic included in the computation. In general, stars with (Li) > 0.5 
are designated secure M31 red giants and stars with (Li) < —0.5 are designated secure MW dwarfs, while stars with 
< (Li) < 0.5 are designated mar ginal M31 red giants and stars with —0.5 < (Li) < are designated marginal MW 
dwarfs (§3.5 of lGilbert et"alll200^ . 

Figure [3] shows the radial velocity distributions of several combinations of stars: only secure M31 RGB stars, secure 
and marginal M31 RGB stars, and secure M31 RGB stars plus all stars that are classified as marginal (M31 red giants 
and MW dwarfs), chosen by use of the diagnostic method both with and without the inclusion of the radial velocity 
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diagnostic. Radial velocity is the most powerful single diagnostic, with a large range of P g i an t /-Pdwarf values (e.g., 
Fig. 10 and §3.4 of lGilbert et al.ll2006l ) and the ability to significantly boost the combined likelihood values. Thus, 
when the radial velocity diagnostic is not included in the overall likelihood calculation, the (Li)y values are in general 
smaller, as can be seen in Figure [16j This causes the significant decrease in the number of secure M31 RGB stars seen 
in the top panel of Figure [3l 

Stars near and more negative than the systemic velocity of M31 have a particularly high probability of being M31 
RGB stars in the radial velocity diagnostic: stars with ^hci = —250 km s _1 have an L v = 2 while stars with i>h e i 
< —300 km s _1 have L v — 5. Consequently, the velocity diagnostic has a large affect on the (Li) values of these 
stars, and they are statistically the most affected by removing the radial velocity diagnostic from the overall likelihood 
calculation. Stars that have velocities in the range over which the radial velocity distributions of the M31 RGB and 
MW dwarf stars overlap (« —200 < «h c i < —125 km s _1 ) will have likelihood values near in the radial velocity 
diagnostic. This means two things: (1) since the L v value is near zero, it adds no power to the (Li) determination, 
and thus the number of objects classifed as secure M31 RGB stars should be similar regardless of whether or not the 
radial velocity diagnostic is used, and (2) since these stars do not have the power of the radial velocity diagnostic, 
they are statistically more likely to land in the marginal (—0.5 <(£,•)< 0.5) regime. The first effect can be seen in 
both the top and middle panels of Figure [31 and the second effect can be seen by comparing the numbers of stars in 
the shaded/dotted histogram in the range —200 < v\ lc \ < —125 km s _1 in the top, middle, and bottom panels. 

The shifting of the overall likelihood distributions to smaller absolute (Li) values due to the removal of the radial 
velocity diagnostic also increases the numbers of marginal MW dwarf stars (and correspondingly decreases the number 
of secure MW dwarf stars). As in the case of the secure M31 RGB stars, removing the radial velocity diagnostic from 
the overall likelihood calculation will have the largest effect on stars that have radial velocities close to or more positive 
than the peak of the dwarf distribution (observed to be at Vh^~ —50 km s -1 in our data set: e.g.. Fig. 2 of Gilbert 
et al. 120061 ). This causes stars that otherwise would be classified as secure MW dwarf stars to fall into the marginal 
(—0.5 < < 0.5) or even secure ((Lj)y> 0.5) M31 RGB regime, depending on the likelihood values of the other 

diagnostics. This boosts the number of stars in this velocity range in the M31 RGB samples selected without the 
inclusion of the radial velocity diagnostic. 

A. 2 Selection of the MSI RGB Sample 

Figure [3] shows that there is no difference in the velocity distributions of secure ((Li)> 0.5) and secure plus marginal 
((Li)> 0.0) M31 RGB stars with t>h c i < —300 km s~ x , chosen with the use of the radial velocity diagnostic. Since 
there is minimal MW dwarf contamination at velocities this negative, we assume that the number of M31 RGB stars 
with (Li)> (465 stars) represents the true number of M31 RGB stars observed in this velocity range. This allows a 
calculation of the level of incompleteness in other samples. There are 446 stars with Vh c \ < —300 km s _1 and (L^)y> 0; 
this sample is 96% complete. Samples selected using the thresholds (Lj)y>[0.5, 0.75, and 1] are [78%, 57%, and 26%] 
complete. 

The vast majority of the dwarf contaminants will have t>h c i > —150 km s _1 . Counts of stars in this velocity range 
for a given sample, compared to an estimate of the expected number of M31 RGB stars in this velocity range, gives an 
estimate of the level of dwarf contamination. The thresholds (i,)y> 0.75 and (Li)y> 1 are quite severe, requiring that 
a star have a probability of being an M31 RGB star that is 5 or 10 times higher than its probability of being an MW 
dwarf star, respectively. These samples are assumed to have minimal dwarf contamination. Based on the numbers of 
stars with Vhd > —150 km s _1 that pass these thresholds (41 stars for (£j)y> 0.75 and 22 stars for (£j)y> 1) and 
the completeness factors calculated above, a complete sample of M31 RGB stars would have approximately 76 stars at 
Vhd > —150 km s _1 . A sample selected using the threshold (Lj)y> 0.5 has 66 stars, which is 87% complete. However, 
the (ij)y > 0.5 sample is expected to be only 78% complete, implying that 9% of the stars with Uhci > —150 km s _1 in 
this sample are MW dwarf star contaminants (overall, this sample has a 1.9% contamination rate). A similar calculation 
for stars selected using the (£j)/> threshold yields a contamination rate of ~ 40% for stars with i>h c i > —150 km s _1 . 
This large local contamination rate corresponds to only a 5.0% contamination rate for the entire sample. A parallel 
calculation using the number of stars with Vh e i > —200 km s _1 yields similar overall contamination rates, confirming 
that the majority of the MW dwarf contamination is at i>h c i > —150 km s _1 . The majority of the dwarf contamination 
comes from blu e stars [panel (d) of Fig. IIP] . since the Nai, CMD, and [Fe/H] diagnostics are all less sensitive for blue 
stars (§4.1.2 of lGilbert etldl^Offfil ). 

The focus of this paper is the substructure discovered along the southeastern minor axis of M31. Therefore, we have 
chosen to use the most complete M31 RGB sample chosen without the use of radial velocity, (Lj)y> 0, as it gives 
the most robust statistics for the cold component, even though it suffers from significant MW dwarf contamination 
at velocities near zero. The removal of the radial velocity diagnostic yields a sample composed of an underlying, 
kinematically unbiased M31 RGB population, as well as a contaminating MW dwarf population at velocities near 
zero. For the analysis presented in this paper, this is preferable to a clean, but kinematically biased M31 RGB sample, 
as the effect of the MW dwarf contamination on the measured properties of the M31 RGB sample can be quantified 
by analysing samples chosen using stricter (Li)y thresholds. 

A. 3 Bias in Measured Spheroid Velocity 
The main effect of the MW dwarf contamination in our chosen M31 RGB sample ((Lj)y> 0) is a positive shift in 
the measured mean velocity of the spheroidal distribution. This shift will affect the measured (i>) sph in the double- 
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Gaussian fit to the combined sample (§[47TJ) , and the (u) sph values from the single-Gaussian fits to individual fields 
(§ !4.2p . The substructure component is cold enough that its measured properties are not significantly affected by the 
dwarf contamination. As discussed above (§A.2), the amount of MW dwarf contamination is small in the (Lj)y> 0.5 
sample and minimal in the (Lj)y> 0.75 sample. Although these samples suffer from significant incompleteness (22% and 
43%, respectively), they are large enough to allow a measurement of the underlying spheroidal component. Maximum- 
likelihood double-Gaussian fits to the (£;}y> 0.5 and (£j)y> 0.75 samples yield (v) sph values of — 302.5lg'g and 
-309.6±s° 7 8 km s"\ respectively. The best-fit a^ h values for the two samples are 123.3±*% and 128.3±g° 7 8 km s -1 , 
respectively. Thus, the best-fit (v) sph values from the (Zi)y> sample should be adjusted by ~ 15 to 20 km s _1 to 
account for the MW dwarf contamination in the sample. This makes the (t>) sph value from the best-fit double Gaussian 
((i;) sph = —287.2 km s _1 ) consistent with the systemic velocity of M31 {v sys — —300 km s^ 1 ). Although the overall 
level of dwarf contamination is expected to increase with radius due to the decreasing surface density of M31 stars, 
there is little difference in the mean offset of (v) sph if the sample is split by radius into inner and outer bins. 
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